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Qh. Abstract 

<D . 

_C • The complete list of definitions for quantities relevant in the anal- 

ysis of SLD/LEP-1 results around the Z-resonance is given. The 
common set of conventions adopted by the programs TOPAZO and 
ZFITTER, following the recommendations of the LEP electroweak 

J3 . working group, is reviewed. The relevance of precision calculations 

is discussed in detail both for pseudo-observables (PO) and for re- 
alistic observables (RO). The model- independent approach is also 
discussed. A critical assessment is given of the comparison between 
TOPAZO and ZFITTER. 
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1 Motivations for the Upgrading of Precision 
Calculations 

The main motivation for upgrading precision calculations around the Z-reso- 
nance with the programs TOPAZO and ZFITTER [|| and for making public 
the results is a reflection of questions frequently asked by the experimental 
community: 

A complete definition of lineshape and asymmetry pseudo observables (POs), 
together with the residual Standard Model (SM) dependence in model-independent 
fits, is needed. This includes a description on what is actually taken from the 
SM. 

Both codes calculate POs. A definition of these POs is needed, showing that 
TDPAZO W and ZFITTER M? use the same definition so that any discrepancy 
is really a measure of missing higher-order corrections. This should include 
quantities like M z , T z , Tf, a® and A® g , and also g ' , g^ A and sin 2 9^ . 

1.1 Goals of this Report 

In 1989 the CERN Report l Z Physics at LEP1', |§ has provided a central doc- 
umentation of the theoretical basis for the physics analysis of the LEP results. 
Although being quite comprehensive, an update on the discussion of radiative 
corrections became necessary in 1995, detailed in the CERN report 'Reports of 
the Working Group on Precision Calculations for the Z Resonance', Q. The 
structure of the latter report was determined by a central part describing the 
situation for the electroweak observables as obtained by various independent 
calculations, including the remaining theoretical uncertainties, followed by com- 
prehensive descriptions of the QCD aspects of electroweak Z physics. 

A new step was taken in early 1998 with a note on the 'Upgrading of Pre- 
cision Calculations for Electroweak Observables' ||, where one focused on the 
calculation of the pseudo-observables. 

It is now time to move a step forward and fully revise the comparisons not 
only for POs but also for realistic observables (ROs), i.e., total cross-sections 
and forward-backward asymmetries, both extrapolated and with realistic cuts. 
Our goal, therefore, has been to upgrade and to compare critically the complete 
TDPAZO and ZFITTER predictions with a particular emphasis on demanding the 
following criteria: 

• Comparisons, after the upgrading, should be consistently better than what 
they were in earlier studies. 

• At the peak all relative deviations among total cross-sections and absolute 
deviations among asymmetries should be below 0.1 per-mill. 



• 



At the wings, typically y^i = M z ± 1.8 GeV, they should be below 0.3 
per-mill. 



It is important to observe that a comparison for ROs at the level of 10~ 4 has 
never been attempted before. 

The numerical results reported in this article are calculated with TOPAZO 
version 4.4 and ZFITTER version 5.20 M. After a careful examination of 
the new upgrading of TOPAZO and ZFITTER contained in these versions we are 
able to report in general a good agreement in our comparisons. The worst case 
for pseudo-observables is represented by the 66-channel. This fact, however, 
was largely expected: this particular channel is where the next-to-leading two- 
loop electroweak corrections are missing and, therefore, here is where we face 
a larger level of theoretical uncertainties. We find satisfactory agreement in all 
the comparisons performed for realistic observables but one: the inclusion of 
initial- final QED interference in the presence of realistic cuts, i.e., acollinearity 
and polar angle cuts and energy thresholds. This fact will be discussed in detail 
in Section 0. 

In this context we would like to emphasise that the implementation of the 
next-to-leading corrections, as well as of any higher-order corrections, makes 
stable all theoretical predictions: the degree of arbitrariness of the various im- 
plementations is reduced with the introduction of newly computed terms. 

Coming back to the reason for the present upgrading, we may say that there 
are additional motivations for it, which we illustrate in the following section. 

1.2 List of Improved I/O Parameters 

For all results, if not stated otherwise we use 

M z = 91.1867 GeV. (1) 

In fixing the set of input parameters we take the lepton masses as in PDG'98 M . 
They are as follows: 

m e = 0.51099907 MeV, m M = 105.658389 MeV, m T = 1.77705 GeV. (2) 

Since all renormalization schemes use the Fermi constant, G F , we refer to a re- 
cent calculation ||, giving an improved value of G F = 1.16637(1) x 10~ 5 GeV - " 
to be compared with the old one, G F = 1.16639(2) x 10~ 5 GeV" 2 j§. 

An important issue concerns the evaluation of oqed at the mass of the Z. 
There is an agreement in our community on using the following strategy. Define 

a{M.) = ^ ss , (3) 

1 Aa(5)(M z ) - A top (M z ) - A t 7 P s (M z ) 

where one has Aa^(M z ) — Aaicpt + Ao^ad- In both codes the input parameter 
is now Aa^, as it is the contribution with the largest uncertainty, while the 
calculation of the top contributions to Aa is left for the code. This should 
become common to all codes. 

The programs TOPAZO and ZFITTER include the recently computed O (a 3 ) 

terms of pi] for Aai op t, and use as default Aa^ = 0.0280398 taken from 



[O. As explained above the latter parameter can be reset by the user. Using 
the default one obtains \/a^\M z ) = 128.877, to which one must add the 
it contribution and the O (aa s ) correction induced by the it loop with gluon 
exchange Q. 

For the improved calculation of Aai Gpt we find the result as reported in 
Tab. (II]). For the t contribution we derive the results listed in Tab.(0). For the 
O (aa s ) corrections induced by the it loop we report the results in Tab.(0) with 
the value of a s (m t ) as reported in Tab.(§). 

2 Analysis of the Measurements 

In the following we take the LEP-1 measurements of hadronic and leptonic cross 
sections and leptonic forward-backward asymmetries as an example to discuss 
the data analysis strategy. 

2.1 The Experimental Strategy 

Technically, each LEP experiment extracts POs, namely M z ,T z ,a1, R e ,fi,T and 



j\0,e,)j,,T ^ gee Section 3 1 f or a definition), from their measured cross-sections and 
asymmetries (realistic observables). The four sets of POs are combined, taking 
correlated errors between the LEP experiments into account, in order to obtain 
a LEP-average set of POs |13| . The latter is then interpreted, for example within 
the frame-work of the Minimal Standard Model. 

Ideally, one would like to combine the results of the LEP experiments at the 
level of the measured cross-sections and asymmetries - a goal that has never 
been achieved so far because of the intrinsic complexity, given the large number 
of measurements with different cuts and the complicated structure of the experi- 
mental covariance matrices relating their errors. As a consequence, the practical 
attitude of the four LEP experiments is to stay with a Model- Independent (MI) 
fit, i.e., from ROs — > POs (©a Standard Model remnant) for each experiment, 
and to average the four sets of POs. The result of this procedure is a set of 
best values for POs which are of course important quantities in their own right. 
The extraction of Lagrangian parameters, M z ,mt, M H , a s (M z ) and a(M z ), is 
based on the LEP-averaged POs. 

There remain several questions to be answered: the main one is, to what 
extend are the POs a Model- Independent (MI) description of the measurements? 
Furthermore, are they valid even in the case where the Standard Model is not the 
correct theory? Note, that many effects are absorbed into the POs. Since POs 
are determined by fitting realistic observables (ROs), one has to clarify what is 
actually taken from the SM (such as imaginary parts and parts which have been 
moved to Z — 7 interference terms and photon-exchange terms) making the MI 
results dependent on the SM. 

The POs are unsatisfactory for many reasons but to some level of accuracy 
they describe well the experimental measurements at the Z peak. How well and 





10 4 x Aaiopt 


TOPAZO 


314.97644 


ZFITTER 


314.97637 



Table 1: 10 4 x A lopt (M z 



m t [GcV] 


168.8 


173.8 


178.8 


TOPAZO 


-0.622230 


-0.585844 


-0.552589 


ZFITTER 


-0.622230 


-0.585844 


-0.552589 



Table 2: 10 4 x A top (M z ). 



a s (M z )/m t [GcV] 


168.8 


173.8 


178.8 


0.116 


-0.108440 


-0.101593 


-0.095371 




-0.108440 


-0.101593 


-0.095371 


0.119 


-0.110994 


-0.103976 


-0.097599 




-0.110994 


-0.103962 


-0.097600 


0.122 


-0.113536 


-0.106347 


-0.099816 




-0.113536 


-0.106347 


-0.099816 



Table 3: TOPAZO (first row) / ZFITTER (second row) results for 10 4 x A t s (M z ) 



a s (M z )/m t [GcV] 


168.8 


173.8 


178.8 


0.116 


0.10631 


0.10589 


0.10548 




0.10631 


0.10589 


0.10548 


0.119 


0.10881 


0.10837 


0.10795 




0.10881 


0.10837 


0.10795 


0.122 


0.11130 


0.11084 


0.11040 




0.11130 


0.11084 


0.11040 



Table 4: TOPAZO (first row) / ZFITTER (second row) results for a s (m t ). 



what is lost in the reduction ROs — ► POs is exactly the kind of question that 
the LEP community is trying to answer. 

In the case of the Standard Model and the measurements of hadronic and 
leptonic cross sections and leptonic forward-backward asymmetries at LEP-1, 
it has been tested by each LEP experiment how the results on the SM pa- 
rameters differ between a SM fit to its own measured ROs, and a SM fit to 
the POs which themselves are derived in an MI fit to the same ROs 113] . In 
the MI fit to determine the POs, the SM initialisation has been performed 
with M z = 91.1867 GeV, m t = 175 GeV, M H = 150 GeV, a s {M 2 z ) = 0.119, 

and Aa£2j(Mf) = 0.02804 {l/a^{M 2 z ) = 128.878). In the two SM fits (to 
ROs and to POs) to be compared, in both cases mt — 173.8 ± 5.0 GeV and 
Aa had( M z) = 0.02804 ± 0.00065 (l/a^(Mf) = 128.878 ± 0.090) are included 
as external constraints. 

For each experiment, the largest difference in central values, relative to the 
fitted errors, is observed for the value of the SM parameter M z , up to 30% 
of the fit error. The fit error itself is unchanged. The shift in central value 
also depends on the Higgs mass used in the SM initialisation of the MI fit to 
determine the POs. One reason for these observations is the following, but more 
detailed studies are needed. The experimentally preferred value of the Higgs- 
boson mass may be different from the value of the Higgs-boson mass used in the 
SM initialisation of the MI fit. In particular, the latter affects the value of the 
Z — 7 interference term for the hadronic cross-section which must be taken from 
the SM for the MI fit. Since the Z — 7 interference term for the hadronic cross- 
section is highly anti-correlated with M z when derived from the LEP-1 data, the 
choice of M H affects the fitted value of M z in the MI fit, and subsequently the 
extraction of SM parameters. For the other four SM parameters, the observed 
differences in fitted central values and errors are usually below 10% to 15% of 
the fit error on this parameter. 

As the LEP average is a factor of two more precise, care has to be taken in 
the averaging procedure. An alternative way to extract SM parameters from the 
LEP measurements, avoiding the intermediate step of POs, is to average directly 
the SM parameters which have been obtained by each experiment through a 
SM fit to its own ROs. This alternative should indeed also be pursued by the 
experiments. 

2.2 The Theoretical Strategy 

Within the context of the SM the ROs are described in terms of some set of 
amplitudes 

A gM = A 7 + A z + non-factorizable, (4) 

where the last term is due to all those contributions that do not factorize into 
the Born-like amplitude, e.g., weak boxes. Once the matrix element A SM is 
computed, squared and integrated to obtain the cross section, a convolution 



with initial- and final-state QED and final-state QCD radiation follows: 

cr(s) = dz H iD (z,s)H &n (z,s)a(zs), (5) 



where H- ln {z 1 s) and Hf[ n (z, s) are so-called, radiator or flux functions accounting 
for Initial- (Final-) State Radiation, ISR (FSR), respectively, and a(zs) is the 
kernel cross-section of the hard process, evaluated at the reduced squared centre- 
of-mass energy s' = zs. 

It is a well-known fact that the structure of the matrix clement changes after 
inclusion of higher-order elcctroweak corrections. One needs the introduction of 
complex-valued form-factors which depend on the two Mandelstam variables s 
and t. The separation into insertions for the 7 exchange and for the Z exchange 
is lost. 

The weak boxes are present as non-resonating insertions to the electroweak 
form- factors. At the Z resonance, the one- loop weak box terms are small, with 
relative contributions < 10~ 4 . If we neglect them, the ^-dependence is turned 
off. The ^-dependence would also spoil factorisation of the form-factors into 
products of effective vector and axial- vector couplings. In all comparisons of 
ROs presented in this report, the weak boxes are taken into account because we 
go off resonance up to y/s = M z ± 3 GeV. 

Full factorisation is re-established by neglecting various terms that are of the 
order O (aT z /M z ). The resulting effective vector and axial- vector couplings are 
complex valued and dependent on s. The factorisation is the result of a variety 
of approximations which are valid at the Z resonance to the accuracy needed, 
and which are indispensable in order to relate POs to ROs. 

After the above mentioned series of approximations we arrive at the so called 
Z-boson pole approximation, which is actually equivalent to setting s = M 2 in 
the form-factors. After de-convoluting ROs of QED and QCD radiation the set 
of approximations transform ROs into POs. 

A source of potential ambiguity is linked to the adopted strategy for extract- 
ing MI POs from the measured ROs. Indeed, the full SM calculation in a MI 
analysis is performed only once at the beginning where one needs to specify in 
addition to M z , which is also a PO, the (remaining) relevant SM parameters 
rrit,M H ,a(M z ),a s (M z ) for the SM-complement of the MI parameterisation, 
RO = RO(PO SM). Subsequent steps in the MI calculation then go directly 
via M z , total and partial widths, and Zff couplings. 

This part of the procedure is particularly cumbersome. However, one has to 
live with the fact that - for practical reasons - the POs will be combined among 
the LEP experiments and survive forever. The cross-section and asymmetry 
measurements will be published by the experiments, but most likely no one will 
ever undertake the effort to combine them. Therefore one is left with the task 
of making sure that the adopted procedure is acceptable. 



3 Pseudo-Observables 

There remains to be investigated the systematic errors arising from theory and 
possible ambiguities in the definition of the MI fit parameters, the POs. 

3.1 Definition of Pseudo-Observables 

Independent of the particular realization of the effective couplings they are 
complex-valued functions, due to the imaginary parts of the diagrams. In the 
past this fact had some relevance only for realistic observables while for pseudo- 
observables they were conventionally defined to include only real parts. This 
convention has changed lately with the introduction of next-to-leading correc- 
tions: imaginary parts, although not next-to-leading in a strict sense, are size- 
able two-loop effects. These are enhanced by factors 7r 2 and sometimes also 
by a factor Nf , with Nf being the total number of fermions (flavour g) colour) 
in the SM. Once we include the best of the two-loop terms then imaginary 
parts should also come in. The latest versions of TDPAZO and ZFITTER therefore 
include imaginary parts of the Z-resonance form factors. 

The explicit formulae for the Zff vertex are always written starting from a 
Born- like form of a pre- factor x fermionic current, where the Born parameters 
are promoted to effective, scale-dependent parameters, 



P f z ln 



[if +ia L ) 7+ -2Q fK is 2 + ia Q }= 7m {Q{ + Q{ 75) , (6) 



where 7+ = 1 + 75 and a Q>L are the SM imaginary parts. Note that imaginary 
parts are always factorized in ZFITTER and added linearly in TDPAZO. 

By definition, the total and partial widths of the Z boson include all cor- 
rections, also QED and QCD corrections. The partial decay width is therefore 
described by the following expression: 

r> = r (z -, fj) = 4c/ r (\G f v \ 2 R f v + \G{\ 2 R{) + A EW/QCD , (7) 

where Cf = 1 or 3 for leptons or quarks (/ = l,q), and the radiator factors 
R v and R A describe the final state QED and QCD corrections and take into 
account the fermion mass m/. 

There is a large body of contributions to the radiator factors in particular for 
the decay Z — > qq; both TOPAZO and ZFITTER implement the results that have 
been either derived or, in few cases, confirmed in some more general setting by 
the Karlsruhe group, see for instance [ ]15| . The splitting between radiators and 
effective couplings follows well defined recipes that can be found and referred to 
in P, n6[ . In particular our choice has been that top-mass dependent QCD cor- 
rections are to be considered as QCD corrections and included in the radiators 
and not in the effective quark couplings. 

The last term, 

EW/QCD x EW/QCD ^ EW ' V u 7 
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accounts for the non-factorizable corrections. The standard partial width, To, 

is 

G M 3 

r ° = 2i7fi = 82 - 945(7) MeV - (9) 

The hadronic and leptonic pole cross-sections are defined by 

r e v l 



r r 

127T ±4\ 

M 2 T 2 

z z 



12tt 



MlTl ' 



(10) 



where T z is the total decay width of the Z boson, i.e, the sum of all partial 
decay widths. Note that the mass and total width of the Z boson are defined 
based on a propagator term \ with an s-dependent width: 

X -\s) = S -Ml+isT z /M z . (11) 

The effective electroweak mixing angles {effective sinuses) are always defined by 

7/ 



4|Q / |sin 2 ^ ff = l-^ 
Re Q 



f 



1 "v 

f 
9a 



where we define 



gf v =Regf v , g{=Reg{. 

The forward-backward asymmetry A is defined via 



(12) 
(13) 



^ FB = 



(14) 



where <r F and a B are the cross sections for forward and backward scattering, 
respectively. Before analysing the forward-backward asymmetries we have to 
describe the inclusion of imaginary parts. A is calculated as 



A =^ 

FB 4 a^, 



(15) 



where 



G F Ml 
V2 



y/PePf QeQ/Re 



a*{Ml)QlQ{ X {s) 



+ 



G\Ml 



p e pMQ e v (ai)lRefe (Q{)*]s\x(s)\ 



(16) 



In case of quark-pair production, an additional radiator factor multiplies <7 VA , 
see also Eq.(p3[). 

This result is valid in the realization where pf is a real quantity, i.e., the 
imaginary parts are not re-summed in pf. In this case 



Qi = Re (gf) + i im (g{) = g { + 1 im (gf) 



g{=if + am(g{). (17) 
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Otherwise Q^ = If ' is a real quantity but pf is complex valued and Eq.fllq 
has to be changed accordingly, i.e., we introduce 



-,/ 



with 



Pf v f > 
(3) 



•,/ 



VpJI} 



(3) 



'/ =-T ~ 2Qy sin 2 0^ . 



(18) 



(19) 



For the peak asymmetry, the presence of p's is irrelevant since they will cancel 
in the ratio. We have 



A of 

FB 



A f = 



A P A 



/. 



2Be[G( r (G{ 



(20) 



The question is what to do with imaginary parts in Eq.(EO). For partial widths, 
as they absorb all corrections, the convention is to use 



I^L I 2 = ( R <J 2 + ( In <. 



(21) 



On the contrary, the PO peak asymmetry A® B will be defined by an analogy of 
equation Eq.(|20|) where conventionally imaginary parts are not included 



A Of 



At 



— A e Af , 



w 



[9 



f) 



(22) 



We note, that Eq. (|22|) is not an approximation of Eq.([20|). Both are POs 
and both could be used as the definition. Numerically, they give very simi- 
lar results: ZFITTER calculates for the two definitions in Eq.([20|) and Eq. (|22|) , 
i™ = 0.0160692 and Af B = 0.0160739. The absolute difference, 0.0000047, is 
more than two orders of magnitude smaller than the current experimental error 
of 0.00096 || . 

In contrast to POs, which are defined, it is impossible to avoid imaginary 
parts for ROs without spoiling the comparison between the theoretical predic- 
tion and the experimental measurement. Then one has to start with Eq.(Hq). 
We will develop Eq.(16) in the realization where imaginary parts are added 
linearly. For the ZZ part of the VA cross-section one derives: 



Re 



G°(S° 



Re 



ot (sir 



(23) 



This collapses to a familiar expression if the axial-vector coefficients are real, 
however one cannot factorize and simplify the p's especially away from the pole 
because of the Zj component. For the Zj part of the VA cross-section one has 
the following result: 



Re[a*(«) X (*)] Re {Q e A Q{) - Im[a*(«) X (*)] Im (G\G^ 



(24) 
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A definition of the PO heavy quark forward-backward asymmetry parameter 
which would include mass effects is 

2g b g b 

Ab = \{*-P)tfJ+Ptf A Y /3 ' (25) 

where (3 is the 6-quark velocity. For A^ TOPAZO/ZFITTER find (with m b = 4.7 
GeV) 

0.102611/0.102634 for m b ^ 0, 

0.102594/0.102617 for m b = 0, (26) 

with a 0.000017/0.000017 difference to be compared with the experimental error 
of 0.0021 [[t3|. The difference is very small, due to an accidental cancellation of 
the mass corrections between the numerator and denominator of Eq. ( |25| ) . This 
occurs for down-type quarks where (g h v ) 2 « (s^ ) 2 /2 and where 



2g e v g e A 2 9 b v 9 b A 

'"■ ~ 4(. 9 -)2 + (^) 2 (g b v ) 2 + (g A ) 2 



A(9 b A ) 2 ll-(9 v ? 



(9 b v ) 2 + (9 b A ) 2 



(27) 



For the c-quark this difference is even bigger (~ 0.000025 for m c = 1.5 GeV, to 
be compared with the experimental error of 0.0044 JL3|), one more example that 
for 6-quarks we meet an accidental cancellation. Note that the mass effect should 
be even smaller since running quark masses seem to be the relevant quantities 
instead of the pole ones. Therefore, our definition of the POs forward-backward 
asymmetry and coupling parameter will be as in Eq.(|22j). 

The most important upgradings in the SM calculation of POs have been 
already described in ||. In particular they consist of the inclusion of higher- 
order QCD corrections, mixed electroweak-QCD corrections fL7| , and next-to- 
leading two-loop corrections of O (a 2 mf) |18| . 

In Rcf. iTq] the two-loop O (a 2 m 2 } corrections are incorporated in the theo- 
retical calculation of M w and sin 2 6*£ . More recently the complete calculation 
of the decay rate of the Z has been made available to us |19] . The only case 
that is not covered is the one of final 6-quarks, because it involves non-universal 
O (a 2 m 2 ) vertex corrections. 

Another development in the computation of radiative corrections to the 
hadronic decay of the Z is contained in two papers, which together provide 
complete corrections of 0(aa s ) to T(Z — > qq) with q — u,d,s,c and b. In 
the first reference of |0 the decay into light quarks is treated. In the second 
one the remaining diagrams contributing to the decay into bottom quarks are 
considered and thus the mixed two-loop corrections are complete. 
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3.2 Model Independent Calculations 

To summarise the MI ansatz, one starts with the SM, which introduces complex- 
valued couplings, calculated to some order in perturbation theory. Next we 
define g^ , g^ as the real parts of the effective couplings and Tf as the physical 
partial width absorbing all radiative corrections including the imaginary parts 
of the couplings and fermion mass effects. Furthermore, we introduce the ratios 
of partial widths 

r t 

(28) 



7? - T " 

H q - — 



r 

Ri = — 



for quarks and leptons, respectively. 

The LEP collaborations report POs for the following sets: 



(M z ,T z ,alR f ,A°J); (M z ,T z ,T h ,gl,g{); (Af z , T z , T, , sin 2 e( s , Pf ). 

(29) 
In order to extract g$ , g^ from Tf one has to get the SM- remnant from Eq. (0) , 
all else is trivial. However, the parameter transformation cannot be completely 
MI, due to the residual SM dependence appearing inside Eq.(0). 

In conclusion, the flow of the MI calculation requested by the experimental 
collaborations is: 



1. Pick the Lagrangian parameters m t ,M H etc. for the explicit calculation 
of the residual SM-dependent part. 



2. Perform the SM initialisation of everything, such as imaginary parts etc. 
giving, among other things, the complement SM. 

3. Select g f v , g{ . 

4. Perform a SM-like calculation of Tf, but using arbitrary values for g' , gf 
and only the rest, namely 



R f R f 



A 



EW/QCD ' 



Im Q-' , Im Q-' 



(30) 



from the SM. 

An example of the parameter transformations is the following: starting from 
M z ,r z ,al,R e ,p,r and A^> T we first obtain 



MS, 



T = MSr 



nl/2 



12nR e 



<Re 



-,1/2 



12tt 



(31) 



With 



Ae ~ V3 



AQ, e 



and 7 



G F M\ 

6V2 7T 



(32) 
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we subtract QED radiation, 



r° 



3 a(M, 



,f2N 



(33) 



and get 



sin 



'efi 



J> 



£2 

7 



v / T~4f-l > 



A 



2sin 2 c e ff 



With 



we further obtain 



2 Q f 



A f = - 



1 

4 + 

4 AH-/ 



(Im^ 



Im q\ 



O J\o 



(34) 



(35) 



sin z 9. 



4 Q/ 



/ 



A f 



~-2\Q f \ S m 2 6 f eS 



(imgtf + (img{y 



(36) 



for / = /j,, r. For quarks one should remember to subtract first non-factorizable 
terms and then to distinguish between R v and R A . 

3.3 Results for Pseudo-Observables 

Having established a common input parameter set (IPS) we now turn to dis- 
cussing the results for pseudo-observables (POs). For POs we use two reference 
sets of values: 

M z = 91.1865 GeV, m t = 171.1 GeV, M H = 76 GeV, a s (M z ) = 0.119, 

(37) 
which corresponds to the minimum of the x 2 of the summer-1998 fit to the SM, 
see [|13|, and 

M z = 91.1867 GeV, m t = 173.8 GeV, M H = 100 GeV, a s {M z ) = 0.119, 

which is our preferred setup in this report. 

For the usual list of POs that enter any SM fit we derive the results of 
Tabs.(J5|-|(]). Here we take into account the updated value for G F and compare 
results calculated with the old value, G F — 1.16639 x 10 -5 GeV -2 , and with 
the new value, G F = 1.16637 x 10 -5 GeV -2 . For all other results presented in 
this report, the new value of G F is used. 
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Observable 


Summer 1998 


01dG F 


New G F 


Diff. 


1/a^ 


128.878 
±0.090 


128.877 




128.877 


M z [GeV] (Input) 


91.1865 


91.1865 




m t [GeV] (Input) 


171.1 


171.1 




M H [GcV] (Input) 


76.0 


76.0 




T z [GeV] 


2.4958 
±0.0024 


2.49543 
2.49564 


2.49538 
2.49559 


0.05 MeV 


4 [nb] 


41.473 
±0.058 


41.4743 
41.4759 


41.4743 
41.4759 


- 


Ri 


20.748 
±0.026 


20.7468 
20.7453 


20.7467 
20.7452 


0.0001 


Ki 


0.01613 
±0.00096 


0.0161823 
0.0161686 


0.0161725 
0.0161588 


0.00001 


J\.Q 


0.1467 
±0.0017 


0.146889 
0.146827 


0.146844 
0.146782 


0.0005 


sin 2 ^ 


0.23157 
±0.00018 


0.231539 
0.231547 


0.231544 
0.231552 


-0.00001 


M w [GeV] 


80.37 
±0.09 


80.3722 
80.3724 


80.3718 
80.3721 


0.4 MeV 



Table 5: Table of POs, first entry is TDPAZO, second is ZFITTER. The ex- 
perimental results show the status of summer 1998 [O. Old G F is G F — 
1.16639 x 10- 5 GeV~ 2 , new G F is G F = 1.16637 x 10~ 5 GeV- 2 . Diff is differ- 
ence between old and new G F . 
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Observable 


Summer 1998 


01dG F 


New G F 


Din". 


l/a® 


128.878 
± 0.0021 


128.877 
128.877 




Af z [GcV] (Input) 


91.1865 


91.1865 




m 4 [GeV] (Input) 


171.1 


171.1 




M H [GcV] (Input) 


76.0 


76.0 




Rb 


0.21590 
±0.00076 


0.215913 
0.215897 


0.215913 
0.215898 


- 


R c 


0.1722 
±0.0048 


0.172223 
0.172224 


0.172222 
0.172223 


- 


A o,b 

FB 


0.1028 
±0.0021 


0.102912 
0.102927 


0.102881 
0.102895 


0.00003 


*% 


0.0734 

±0.0045 


0.0735700 
0.0735365 


0.0735456 
0.0735121 


0.00002 


Ab 


0.935 
±0.018 


0.934724 
0.934678 


0.934720 
0.934674 


- 


Ac 


0.668 
±0.028 


0.667806 
0.667784 


0.667787 
0.667765 


0.00002 



Table 6: Table of POs, first entry is T0PAZ0, second is ZFITTER. The ex- 
perimental results show the status of summer 1998 [Of. Old G F is G F — 
1.16639 x 10- 5 GeV~ 2 , new G F is G F = 1.16637 x lO^Ge^- 2 . Diff is differ- 
ence between old and new G F . 
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The full list contains more POs and is given in Tab. (pi), where we include 
the relative and absolute difference between TOPAZO and ZFITTER in units of 
per-mill: 

T — Z 

S = 10 3 x — — A = 10 3 x (T - Z) . (39) 

In Tab.(|7j) we report also some POs which are not usually taken into account in 
fitting the experimental data. With the exception of sin 2 6 h oS , for which we find 
a difference of 0.4 per-mill, the relative deviation is always (well) below 0.15 
per-mill. The larger difference in sin b B g is perhaps not surprising since the b 
sector did not undergo any update aimed to including next-to-leading two-loop 
electroweak effects in mt, which are not available for this channel. 

In Tabs.(||-^|) we analyse the POs as a function of M H in logarithmic spacing, 
M H = 10, 30, 100, 300, 1000 GeV, therefore including the region of low M H (here 
M z = 91.1867 has been used). 

The variations of POs as a function of the Higgs boson mass are an important 
issue related to the theoretical uncertainty in the determination of constraints 
on SM parameters. From the most recent study |ll| one derives that the region 
below a Higgs mass of about 70 GeV has a comparatively larger theoretical 
uncertainty, although the current 90 GeV lower limit on the Higgs mass from 
the direct search makes it less interesting. 

In Tabs.(&-0) we compare a relevant set of POs in the preferred calculational 
setup of TOPAZO and ZFITTER. In Tab.(||) we report the relative deviations, in 
per-mill, between the two predictions. Everywhere this deviation is (well) below 
0.1 per-mill, even at very low values of the Higgs mass. For the asymmetries 
of Tab.(^|) we report absolute deviations in units of 10~ 3 . The largest absolute 
deviation is found for A e at M H = 10 GeV, giving 0.14 x 10~ 3 . For M H = 
100 GeV all absolute deviations in Tab.(||) are below 0.08 x 10~ 3 . 

The good agreement of POs calculated by TOPAZO and ZFITTER verifies a 
posteriori the consistency of definitions for POs in the two programs, although 
one should understand that the agreement is necessary but not sufficient for 
consistency. We come back to this question when we discuss realistic observables 
and their calculations in terms of POs. 

It is instructive to compare few examples with the old results of Q , obtained 
for M H = 300 GeV. For sin 2 Q*g the relative difference T/Z has moved from 0.21 
to 0.004 per-mill at M H = 300 GeV and it is now everywhere below 0.06 per- 
mill, reached at very low Higgs masses. For M w it was 0.087 per-mill and it 
is now 0.015 per-mill with a maximum of 0.029 for very large values of M H . 
Finally, the absolute deviation for A e was 0.38 x 10~ 3 and it is now 0.05 x 10 -3 
with a maximum of 0.14 x 10 -3 at low values of M„. 
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Observable 


TDPAZO 


ZFITTER 


10 3 x T T Z 


l/a( 5 )(M z ) 


128.877 


128.877 




l/a(M z ) 


128.887 


128.887 




M w [GeV] 


80.3731 


80.3738 


-0.009 


o\ [nb] 


41.4761 


41.4777 


-0.04 


a°[nb] 


1.9995 


1.9997 


-0.12 


r fc [Gov] 


1.74211 


1.74223 


-0.07 


Y z [GcV] 


2.49549 


2.49573 


-0.10 


r„ [Mcv] 


167.207 


167.234 


-0.16 


T e [MeV] 


83.983 


83.995 


-0.14 


IV [MeV] 


83.983 


83.995 


-0.14 


T r [MeV] 


83.793 


83.805 


-0.14 


T u [MeV] 


300.129 


300.154 


-0.08 


T d [MeV] 


382.961 


382.996 


-0.09 


T c [MeV] 


300.069 


300.092 


-0.08 


T b [MeV] 


375.997 


375.993 


0.01 


r inv [GcV] 


0.50162 


0.50170 


-0.16 


Ri 


20.7435 


20.7420 


0.07 


K 


0.215829 


0.215811 


0.08 


R° c 


0.172245 


0.172246 


-0.01 


sin 2 ^ 


0.231596 


0.231601 


-0.02 


sin 2 ^ ff 


0.232864 


0.232950 


-0.37 


sin 2 ^ ff 


0.231491 


0.231495 


-0.02 


Pe 


1.00513 


1.00528 


-0.15 


Pb 


0.99413 


0.99424 


-0.11 


Pc 


1.00582 


1.00598 


-0.16 


Observable 


TDPAZO 


ZFITTER 


10 3 x (T - Z) 


^ 


0.016084 


0.016074 


0.01 


A 0,b 

FB 


0.102594 


0.102617 


-0.02 


^ FB C 


0.073324 


0.073300 


0.02 


^e 


0.146440 


0.146396 


0.04 


A b 


0.934654 


0.934607 


0.05 


A c 


0.667609 


0.667595 


0.01 



Table 7: Complete table of POs, from T0PAZ0 and ZFITTER. 
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M H [GeV] 


Observable 


10 


30 


100 


300 


1000 


T z [GeV] 


2.49298 


2.49618 


2.49549 


2.49227 


2.48732 




2.49322 


2.49645 


2.49573 


2.49240 


2.48751 




-0.10 


-0.11 


-0.10 


-0.05 


-0.08 


al [nb] 


41.4739 


41.4744 


41.4761 


41.4788 


41.4830 




41.4748 


41.4761 


41.4777 


41.4798 


41.4831 




-0.02 


-0.04 


-0.04 


-0.02 


-0.002 


<r? [nb] 


1.99797 


1.99851 


1.99947 


2.00062 


2.00209 




1.99811 


1.99874 


1.99970 


2.00074 


2.00208 




-0.07 


-0.12 


-0.12 


-0.06 


-0.005 


Ri 


20.7580 


20.7527 


20.7435 


20.7330 


20.7199 




20.7570 


20.7511 


20.7420 


20.7322 


20.7200 




+0.05 


+0.08 


+0.07 


+0.04 


-0.005 


sin^f 


0.230698 


0.231044 


0.231596 


0.232175 


0.232845 




0.230712 


0.231056 


0.231601 


0.232176 


0.232838 




-0.061 


-0.052 


-0.022 


-0.004 


+0.030 


M w [GeV] 


80.4587 


80.4298 


80.3731 


80.2989 


80.2045 




80.4583 


80.4297 


80.3738 


80.3001 


80.2068 




+0.005 


+0.001 


-0.009 


-0.015 


-0.029 


R b 


0.215759 


0.215794 


0.215829 


0.215839 


0.215824 




0.215763 


0.215775 


0.215811 


0.215845 


0.215857 




-0.02 


+0.09 


+0.08 


-0.03 


-0.15 


R c 


0.172305 


0.172280 


0.172245 


0.172213 


0.172184 




0.172301 


0.172281 


0.172246 


0.172210 


0.172174 




+0.02 


-0.01 


-0.01 


+0.02 


+0.06 



Table 8: Table of POs as a function of the Higgs boson mass; first entry is 
T0PAZ0, second is ZFITTER, third entry is 1-T/Z in pcr-mill. 
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M H [GeV] 


Observable 


10 


30 


100 


300 


1000 


A 0A 

FB 


0.017672 


0.017052 


0.016084 


0.015098 


0.013994 




0.017647 


0.017031 


0.016074 


0.015096 


0.014006 




+0.03 


+0.02 


+0.01 


+0.002 


-0.01 


A 0,b 

FB 


0.107564 


0.105656 


0.102594 


0.099373 


0.095637 




0.107587 


0.105665 


0.102617 


0.099410 


0.095711 




-0.02 


-0.01 


-0.02 


-0.04 


-0.07 


A 0,c 

FB 


0.077216 


0.075714 


0.073324 


0.070827 


0.067949 




0.077157 


0.075663 


0.073300 


0.070824 


0.067983 




+0.06 


+0.05 


+0.02 


+0.003 


-0.03 


•A-e 


0.15350 


0.15078 


0.14644 


0.14188 


0.13660 




0.15340 


0.15069 


0.14640 


0.14187 


0.13666 




+0.10 


+0.09 


+0.04 


+0.01 


-0.06 


A b 


0.935220 


0.935003 


0.934654 


0.934283 


0.933837 




0.935165 


0.934947 


0.934607 


0.934251 


0.933844 




+0.06 


+0.06 


+0.05 


+0.03 


-0.01 


J\c 


0.670709 


0.669517 


0.667609 


0.665601 


0.663267 




0.670666 


0.669481 


0.667595 


0.665605 


0.663302 




+0.04 


+0.04 


+0.01 


-0.004 


-0.04 



Table 9: Table of POs as a function of the Higgs boson mass; first entry is 
T0PAZ0, second is ZFITTER, third is 10 3 X (T - Z). 
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3.4 Theoretical Uncertainties for Pseudo-Observables 

Here we discuss the theoretical uncertainties associated with the SM calculation 
of POs. In Tabs.(J10|-[ll|) we give the central value, the minus error and the plus 
error as predicted by TOPAZO and compare with the current total experimental 
error where available. The procedure is straightforward: both codes have a 
preferred calculational setup and options to be varied, options having to do with 
the remaining theoretical uncertainties and the corresponding implementation 
of higher-order terms. To give an example, we have now LO and NLO two-loop 
EW corrections but we are still missing the NNLO ones and this allows for 
variations in the final recipe for pf, etc. 

TOPAZO has been run over all the options remaining after implementation of 
NLO corrections, and all the results for POs are collected. We use 

- central for the value of the PO evaluated with the preferred setup; 

- minus error for VO central - min opt PO; 

- plus error for max opt PO - PO centr ai- 



Observable 


central 


minus error 


plus error 


total exp. error 


l/a^(M z ) 


128.877 


- 


- 




l/a(M z ) 


128.887 


- 


- 




M w [GeV] 


80.3731 


5.8 MeV 


0.3 MeV 


64 MeV 


<r° [nb] 


41.4761 


1.0 pb 


1.6 pb 


58 pb 


a\ [nb] 


1.9995 


0.17 pb 


0.26 pb 


3.5 pb 


r„ [MeV] 


167.207 


0.017 


0.001 




T e [MeV] 


83.983 


0.010 


0.0005 


0.10* 


T M [MeV] 


83.983 


0.010 


0.0005 




r r [MeV] 


83.793 


0.010 


0.0005 




T u [MeV] 


300.129 


0.047 


0.013 




T d [MeV] 


382.961 


0.054 


0.010 




T c [MeV] 


300.069 


0.047 


0.013 




T b [MeV] 


375.997 


0.208 


0.077 




Thad [GeV] 


1.74211 


0.26 MeV 


0.11 MeV 


2.3 MeV* 


Tinv [GeV] 


0.50162 


0.05 MeV 


0.002 MeV 


1.9 MeV* 


T z [GeV] 


2.49549 


0.34 MeV 


0.11 MeV 


2.4 MeV 



Table 10: Theoretical uncertainties on POs from TOPAZO. The experimental 
error is that of summer 1998 |L3|. *) assumes lepton universality. 
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Observable 


central 


minus error 


plus error 


total exp. error 


Ri 
R° b 
R° c 


20.7435 
0.215829 
0.172245 


0.0020 
0.000100 
0.000005 


0.0013 
0.000031 
0.000024 


0.026* 
0.00074 
0.0044 


sin 2 ^ 
sin 2 ^ ff 
sin 2 ^ ff 


0.231596 
0.232864 
0.231491 


0.000035 
0.000002 
0.000029 


0.000033 
0.000048 
0.000033 


0.00018* 


FB 


0.016084 
0.102594 
0.073324 


0.000057 
0.000184 
0.000142 


0.000060 
0.000195 
0.000149 


0.00096* 
0.0021 
0.0044 


A b 

A c 


0.146440 
0.934654 
0.667609 


0.000259 
0.000032 
0.000114 


0.000275 
0.000005 
0.000103 


0.0017* 
0.035 
0.040 


Pe 
Pb 
Pc 


1.00513 
0.99413 
1.00582 


0.00010 
0.00048 
0.00010 


0.000005 
0.000001 
0.000005 


0.0012* 



Table 11: Theoretical uncertainties on POs from T0PAZ0. The experimental 
error is that of summer 1998 p3|. *) assumes lepton universality. 



One can see from Tabs . (pL 0|-|l 1[) that there is a sizeable reduction of the the- 
oretical uncertainty associated with POs compared to the findings of pj. This 
is in accordance with the work of l| and is mainly due to the implementation 
of next-to-leading corrections in T0PAZ0 and ZFITTER. We do not show any esti- 
mate for theoretical uncertainty in POs from ZFITTER, since they are typically 
more narrow. 

Within T0PAZ0 the central values are defined by the following flags: 0U0= ' S ' 
(fixed), 0U1='Y', 0U2='N' (fixed), 0U3='Y' (fixed), 0U4='N' , 0U5='Y', 0U6='Y', 
0U7= ' N ' , 0U8= ' C ' . Plus and minus errors are obtained by changing the flags to 
the following values: 0U1='N', 0U4='Y', 0U5='N', 0U6='N', 0U7='Y>, 0U8='L' 
or 0U8='R' (one by one). ZFITTER's central values are produced with the default 
flag setting. Plus and minus errors are obtained by changing the flags to the 
following values: SCAL=0,4; HIGS=0,1; SCRE=0,1,2; EXPR=0,1,2; EXPF=0,1,2; 
HIG2=0,1. 

When performing a SM analysis of measured POs, several SM fits should 
be performed, changing the flags as indicated above. The differences in fitted 
values are a measure of the theoretical uncertainty in extracting SM parameters 
from the measured POs. 
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4 Realistic Observables 

The ROs, measured cross-sections and forward-backward asymmetries, are com- 
puted in the context of the SM, see pi . Thus the comparison between TDPAZO 
and ZFITTER. is mainly a SM comparison. In addition, however, one of the goals 
will be to pin down 

• the definition of POs; 

• the calculation of ROs in terms of the defined POs for the purpose of MI 
fits, showing that for POs with values as calculated in the SM, the ROs 
are by construction identical to the full SM RO calculation. 

The last point requires expressing p's and effective mixing angles in terms 
of POs, assuming the validity of the SM. After this transformation the ROs 
will be given as a function of the POs at their SM values. This is not at all a 
trivial affair because of gauge invariance and one should remember that gauge 
invariance at the Z pole (on-shell gauge invariance) is entirely another story 
from gauge invariance at any arbitrary scale (off-shell gauge invariance). Some 
of the re-summations that are allowed at the pole and that heavily influence 
the definition of effective Z couplings are not trivially extendible to the off- 
shell case. Therefore, the expression for RO=RO(PO), at arbitrary s, requires 
a careful examination and should be better understood as RO=RO(PO,SM), 
that is, for example: 

ctmi = o-sm (Rt,A£,--'->gt,g{-> Pf, sin 2 e f cS ; residual Sm) . (40) 

As long as the procedure does not violate gauge invariance and the POs are 
given SM values, there is nothing wrong with the calculations. It is clear that 
in this case the SM ROs coincide with the MI ROs. 

The next question is, of course, do the ROs in MI calculations agree for POs 
not having SM values - at least over a range of PO values corresponding to 
current experimental errors on POs for a single experiment, i.e., two to three 
times the error on the LEP-average POs? It is clear that the present procedure, 
SM fixed and POs varying around their SM values, is wrong in principle but 
one should content oneself with testing its accuracy. 

There is another reason to be worried, one should avoid any interpretative 
strategy such that the pattern becomes: raw data — > RO decoded into PO — * 
Lagrangian parameters (any Lagrangian, SM, MSSM, etc.), if only one decoder 
(code that allows for MI studies) has been built for that purpose. The glimpse 
we want to have of nature should not depend on the decoder. 

4.1 Setup for Realistic Observables 

In this note we discuss ROs for s-channel processes, thus excluding Bhabha 
scattering. The ROs, cross-section and forward-backward asymmetries, are 



24 



INPUT 


M z 


= 91.1867, mt = 173.8, M H = 100,a s (M z ) = 0.119 


\fs 




M z - 3, M z - 1.8, M z , M z + 1.8, M z + 3 


Sq/s 




0.01 


WEAK BOXES 




YES 


IFI 




NO (-» Section g) 


ISPP 




NO (-> Section 0) 



Table 12: Extrapolated setup for the calculations of ROs. Masses and y/s are 
in GeV. 



computed for the setup specified in Tab.(|12|). This setup will be referred to 
as the fully extrapolated setup. 

The choice of the energies is dictated by the fact that the most precise 
SLD/LEP-1 measurements are at the pole and at ±1.8 GeV away from the pole. 
The parameter sq is usually referred to as the s'-cut, i.e., so = s' min , where the 
following definition applies: 
VV is the centre- of-mass energy of the e + e~ system after initial-state radiation. 

In general the effects of initial-final QED interference (IFI) and of initial- 
state pair-production are not included. They will be discussed separately in 
Sections [Z] and [h], respectively. 

When SM parameters are varied we use 



M„ 



10,30,100,300,1000 GeV. 



4.2 Next-to-Leading and Mixed Corrections for Realistic 
Observables 

The inclusion of mixed two-loop correction for RO, at s ^ M z , can only repre- 
sent an approximation to the real answer. Consider the term giving 



EW/QCD 



= N c 



a a crM, 



2sl 



it 12ir 
fl, 



(.9 9 + ) 4 + (g q ) 4 



C P A? ( 



+6/, 



(3) £« 

se 



C'A 
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with g^ = g% ±g q A and C F = 4/3. The functions A^ are given by an expansion in 
a s and, in computing the Z width one sets s = M z obtaining a gauge invariant 

answer. For arbitrary s the following happens.[J A;, (x) is very simple, it quickly 

1 A. Czarnccki and J. Kiihn, private communication. 
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approaches its asymptotic value —3/8 as x grows. So it is legitimate to use the 
same value as at x = 1, i.e., — 0.37± 0.04. 

With A3 ' the situation is a bit more complicated. Its value can of course 
still be found using the formula (14) of the first Ref. of jffl], (although the error 
bar will be larger than at x = 1). However, for a heavy off-shell Z the decay into 
a real W (plus a pair of quarks) becomes more important, and the coefficient 
A 3 ' gives only one part of the full mixed QCD/electroweak corrections. An 
estimate of the QCD corrections to the real W emission is given in [E0[. The 
fact that H^(on-shell)qg is a genuine four-fermion event does not help too much: 
the whole issue of separating two- from four-fcrmion events at sufficiently high 
energy has not yet been systematised. 

In addition, if we are away from the Z pole, we have to take into account 
7-exchange. 

The strategy adopted by T0PAZ0 in this case will be the following: the 
amplitude squared due to Z-exchange is something like 

sf(-s) 



(s- M 2 ) + s 2 T 2 z / 'Ml 
where k denotes a collection of coupling constants. This we rewrite as 



(42) 



sf{s)-Mlf{Ml) T\Ml fjMl) 



(s ~ Ml ) +s 2 T 2 z I Ml (s - Ml ) + s 2 T 2 z J Ml r 2 

The splitting is motivated by the fact that f(s) is not gauge invariant while 
f(M 2 ) is. Moreover, for on-shell Z bosons we have at our disposal an improved 
calculation, e.g., including mixed two- loop effects. Thus we write 

sf(s)-M 2 f(M 2 ) T 2 M 2 

■F imp (M 2 ), (44) 



{s-MlY + s 2 Y\ J M 2 z (s - M 2 z ) +s 2 T 2 z /M 2 z 

with F — nf /T 2 , and i*i mp is the two-loop corrected, on-shell, expression. 

Within ZFITTER the implementation of CKHSS correction was done in a sim- 
plified way. The numbers for non-factorized O {aa s ) corrections for different 
channels Z — v //, reported in |l7|), are hard-wired to the code for calculating 



POs. An analogy of Eq.(44) was used for ROs with the inclusion of, prop- 
erly normalised, non-factorized O (aa s ) corrections. A detailed comparison of 
T0PAZ0/ZFITTER numbers with/without CKHSS corrections (DD versus DDD) 
was done and excellent agreement was found. It does not look surprising since 
the correction is small (~ 0.3 per-mill) and its crude implementation works in 
practice. 

The same applies for the next-to-leading, O (a 2 m 2 ) corrections [[jlj. Here 
TDPAZO uses 

sf(s) _ */ LO ( S )-Af z 2 / LO (Af2) 



K ■ 



Mlf + s 2 T 2 z I M 2 z (s - M 2 Z ) 2 + s 2 T 2 z J M 2 Z 



2 1\,t2 



ViM 



( S -MlY + S 2 T 2 z /M\ 
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■F NLO (Af 2 ), (45) 



with LO,NLO indicating leading and next-to-leading corrections. 

Implementation of NLO, two-loop EW corrections, into ZFITTER is very in- 
volved and cannot be described in all details here. It makes use of the FORTRAN 
code m2tcor [plf and is due to common work with G. Degrassi and P. Gam- 
bino done in February 1998 p2fl . The full collection of the relevant formulae 
is presented in |J6fl . Again, a careful comparison of TOPAZO/ZFITTER numbers 
with/without NLO corrections was undertaken and good agreement was reg- 
istered. Since the implementation into TOPAZO is completely independent, the 
agreement is a convincing argument to conclude that both implementations are 
correct. 

4.3 Final-State Radiation 

One should realize that s' is not equivalent to the invariant mass of the final-state 
// system, M 2 (//), due to final state QED and QCD radiation. Furthermore, 
in the presence of a s'-cut the correction for final state QED radiation is simply 



«S|dW = 7 0/— ■ («) 



while for a M 2 -cut the correction is more complicated, see |23 . For full angular 
acceptance one derives the following corrections: 



a(s) = ^-Q}a Q ( S ){ -2.r- 



z + ^- + 2lnx)\nA T (47) 

2 ) mi 



+z(l + ~) lnz + 2C(2) -2Li 2 (x) -21nir + ^ -Zz-^- 

Here we have introduced 

x = l-z, z = M 2 (//) /s. (48) 

For an s'-cut both QED and QCD final-state radiation are included through an 
inclusive correction factor. For M 2 -cut the result remains perfectly defined for 
leptons, however for hadronic final states there is a problem. This has to do 
with QCD final-state corrections. Indeed we face the following situation: 

• for e + e~ — ► ffj the exact correction factor is known at O (a) even in the 
presence of a M 2 (fJ) cut ||; 

• the complete set of final-state QCD corrections are known up to O (o^ ) 
(see J24| and also []15J) only for the fully inclusive setup, i.e., no cut on the 
// invariant mass; 

• the mixed two-loop QED/QCD final-state corrections are also known only 
for a fully inclusive setup P5| : 

• at O (a s ) QCD final-state corrections in presence of a M 2 -cut follow from 
the analogous QED calculation |26[|. 
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The ideal thing would be to have QED QCD final-state radiator factors 
R, with a kincmatical cut imposed on M 2 (ff). Missing this calculation, that 
would give the correct O (aa s ) factors with cuts, we have three options 

pFS _ 1 | nFS i pFS 

n — 1 "+" -"■QED.cut "T ^QCD.cxt i 

^ = (1 + ^QCD,ext) (1 + -RQED,cut) ; 

R = 1 + i? QED . cut + ^QCD.cut • ( 49 ) 



The (QCD,ext) corrections are understood up to 0(a|), while those corre- 
sponding to the (QCD, cut) setup are only computed at O (a s ). The first of 
Eq.(ff9[) is our preferred option. 

5 De-Convoluted Realistic Observables 

Our goal in describing the theoretical uncertainties for realistic observables is 
twofold. First we want to discuss the effect of QED radiation by comparing dif- 
ferent radiators and then we have to give a critical assessment of the theoretical 
uncertainty in the predicted cross-sections, de-convoluted of QED effects, i.e., 
the purely weak uncertainty. We therefore define several levels of de-convolution: 



• 



• 



• 



Single-de-convolution (SD), giving the kernel cross-sections without initial- 
state QED radiation, but including all final-state correction factors. 

Double-de-convolution (DD), giving the kernel cross-sections without ini- 
tial- and final-state QED radiation and without any final-state QCD ra- 
diation. There is an additional level, to be called DDD, and the difference 
between DD and DDD deserves a word of comment. The improvement 
upon naive electroweak/QCD factorisation contains two effects, the FTJR 
correction |]27[ which gives the leading two-loop answer for the 66-channcl 
and the CKHSS correction JL7J which gives the correct answer for the 
remaining mixed corrections in all quark channels. 

In DD-mode FTJR and CKHSS corrections are kept while in DDD-mode 
they are excluded. This option allows us to keep under control the imple- 
mentation of the new CKHSS correction. 

DD, DDD with only Z - Z exchange (DDZ, DDDZ), weak boxes are not 
included, 

DDD with only Z(Bj (DDZG), i.e. noZ-7 interference and weak boxes 
are not included. 



Rather than comparing only the complete results for ROs, i.e., including 
initial-state QED radiation, final-state QED and QCD radiation, initial-state 
pair-production and initial-final QED interference, we do more. The reason is 
given by the observation that often an agreement on the complete result may be 
a consequence of several compensations of the single components. This is why 
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we want to compare component by component and the procedure will allow us 
to formulate a quantitative statement on the overall theoretical uncertainty. In 
this way the errors are independent of the amount of cancellation between the 
various components. 
We define 

a T (X) and A (X) = °" FB / ; , with X= SD,DD,DDD,DDZ,DDZG. 

<J T (X) 

(50) 
For SD-quantities that contain final-state QED radiation, we further distinguish 
between two series, the so-called M 2 -series where a cut on M 2 (ff) is applied 
and the so-called s'-series. For SD setup the latter implies that no cut is ap- 
plied (fully extrapolated setup) . In this way also the final-state QED correction 
factors with/without cuts on the final-state fermions can be compared. 

5.1 De-Convoluted Cross-Sections 

By comparing SD with DD quantities we are able to disentangle the effect of 
initial-state QED radiation from final-state QED QCD radiation. 

The DDZ or DDZG modes are included for convenience of the reader but 
deserve an additional comment. Clearly, away from the Z-peak, diagrams with 
Z or 7 exchanges are not gauge invariant and, therefore, we expect deviations in 
the result of the two codes. However, they are useful to understand the pattern 
of agreement and also to show the internal consistency of codes, for instance 



°7 \ M z ) = °7 



3 a(M 

4 7T 



2^-1 r 



4 } -k 



(51) 



is an equality between RO and PO that must be satisfied. The comparison 
between DDZ and DDZG modes, moreover, gives an estimate of the Z — 7 
interference effects, before folding with QED radiation. 

In Tab.(|l3|) we start our comparison for RO de-convoluted observables show- 
ing results for <jf in DD and SD (no-cut and M 2 -cut) DDD modes. The reference 
point has been fixed to M H = 100 GeV. The relative deviations 5tz =T/Z-1 
in per-mill are shown in Fig. ffl. At ^/s = M z we register 0.05 per-mill for 
muons both for DD mode and for SD, no-cut mode. For hadrons we have —0.02 
per-mill in DD mode and —0.07 per-mill in SD, no-cut mode. There are tiny 
variations when we consider the SD M 2 -cut mode. 

The differences between DD and SD describe different implementations of 
QED final-state radiation and of QCD corrections. Our comparison between 
the two SD branches shows that final-state QED correction factors are correctly 
implemented for the fully inclusive setup (no-cut) and for a cut on the invariant 
mass of the // pair. The agreement between DD and SD for hadrons shows 
that also final-state QCD factors are under control. 

As we have already illustrated there are other sources of final-state a s - 
dependent corrections, due to the non-factorisation of QCD and purely elec- 
troweak effects. The SD mode also accounts for non-factorizable two-loop ef- 
fects and the agreement between two calculations in SD mode is an agreement 
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for the R VtA factors of Eq.(M) and for the interplay between them and the 
FTJR/CKHSS effects. From Fig. | we observe a -0.02, -0.07 and +0.07 per- 
mill differences for u, d and b quarks at the peak. Therefore the overall agreement 
for the hadronic cross-section, —0.02 per-mill, is also the result of some partial 
compensation between contributions from up- and down- type quarks. 

Note that the agreement remains very good also for 6-quarks where next- 
to-leading corrections are not available and where, therefore, one would expect 
larger deviations between T0PAZ0 and ZFITTER. 

Note the following relation between the pseudo-observable Rb and the ratio 
of DD cross-sections 



Rb — 

Chad 



0.00146 T0PAZ0 , . 

V^=M 1 0.00146 ZFITTER. 



The difference reflects the S'M-remnant effect, since the ratio of RO cross- 
sections has 7-exchange, imaginary parts, . . ., and (substantially negligible) 
weak boxes. 
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Centre-of 


-mass energy in GeV 




M z -3 


M z - 1.8 


M z 


M z + 1.8 


M z +3 


a^ [nb] DD = DDD 


0.29999 


0.65718 


2.00341 


0.65856 


0.31047 




0.30003 


0.65724 


2.00331 


0.65863 


0.31051 


an [nb] SD no- cut 


0.30055 


0.65839 


2.00711 


0.65978 


0.31104 




0.30058 


0.65845 


2.00700 


0.65985 


0.31108 


a^ [nb] SD - M 2 -cut 


0.30047 


0.65821 


2.00656 


0.65960 


0.31095 




0.30052 


0.65832 


2.00659 


0.65971 


0.31102 


a u [nb] DD 


0.99648 


2.21803 


6.82893 


2.23330 


1.04203 




0.99682 


2.21861 


6.82913 


2.23355 


1.04214 


o u [nb] SD - no-cut 


1.04290 


2.32118 


7.14541 


2.33633 


1.08993 




1.04330 


2.32183 


7.14551 


2.33647 


1.08996 


a u [nb] SD - M 2 -cut 


1.04277 


2.32090 


7.14452 


2.33604 


1.08979 




1.04320 


2.32162 


7.14486 


2.33626 


1.08986 


a d [nb] DD 


1.26996 


2.84741 


8.79775 


2.86549 


1.32868 




1.27040 


2.84820 


8.79838 


2.86600 


1.32893 


<Jd [nb] SD - no-cut 


1.31395 


2.94596 


9.10195 


2.96450 


1.37458 




1.31444 


2.94683 


9.10268 


2.96502 


1.37482 


a d [nb] SD - M 2 -cut 


1.31391 


2.94587 


9.10166 


2.96441 


1.37453 




1.31441 


2.94672 


9.10249 


2.96496 


1.37479 


a c [nb] DD 


0.99648 


2.21803 


6.82893 


2.23330 


1.04203 




0.99682 


2.21861 


6.82913 


2.23355 


1.04214 


er c [nb] SD - no-cut 


1.04267 


2.32070 


7.14397 


2.33588 


1.08972 




1.04307 


2.32133 


7.14405 


2.33601 


1.08975 


a c [nb] SD - M 2 -cut 


1.04262 


2.32057 


7.14359 


2.33575 


1.08967 




1.04303 


2.32124 


7.14375 


2.33592 


1.08971 


cr b [nb] DD 


1.25204 


2.80753 


8.67631 


2.82663 


1.31091 




1.25226 


2.80787 


8.67577 


2.82681 


1.31101 


Ob [nb] SD - no-cut 


1.28945 


2.89170 


8.93787 


2.91233 


1.35080 




1.28995 


2.89266 


8.93909 


2.91306 


1.35115 


a h [nb] SD - M 2 -cut 


1.28945 


2.89168 


8.93780 


2.91231 


1.35079 




1.28995 


2.89266 


8.93908 


2.91306 


1.35115 


CThad [nb] DD 


5.78492 


12.93841 


39.92967 


13.02421 


6.05233 




5.78670 


12.94148 


39.93079 


13.02591 


6.05313 


Chad [nb] SD - no-cut 


6.00291 


13.42550 


41.43114 


13.51353 


6.27961 




6.00518 


13.42948 


41.43401 


13.51559 


6.28050 


cr had [nb] SD - M 2 -cut 


6.00265 


13.42490 


41.42929 


13.51293 


6.27933 




6.00500 


13.42907 


41.43272 


13.51517 


6.28030 



Table 13: T0PAZ0/ZFITTER comparison of 07 de-convoluted: DD and SD fully 
extrapolated mode and M 2 -cut of 0.01s. Here M H = 100 GeV. 
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Figure 1: Relative deviations between TDPAZO and ZFITTER. for total fermionic 
cross-section. DD-mode is shown with SD-modes corresponding to no-cut or a 
M 2 (//)-cut of 0.01s. 
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5.2 De-Convoluted Asymmetries 

In Tab.Jl4|) we show DD and SD modes for the de-convoluted muonic and heavy 
quark forward-backward asymmetries. 







Centre-of- 


mass energy in GeV 




M z - 3 


M z - 1.8 


M z 


M z + 1.8 


M z + 3 


A» DD 


-0.26170 


-0.15037 


0.01745 


0.17510 


0.27002 




-0.26167 


-0.15037 


0.01741 


0.17502 


0.26991 


A^ B SD - no-cut 


-0.26122 


-0.15010 


0.01742 


0.17478 


0.26952 




-0.26119 


-0.15010 


0.01737 


0.17469 


0.26941 


A* SD - M 2 -cut 


-0.26128 


-0.15013 


0.01742 


0.17481 


0.26958 




-0.26122 


-0.15011 


0.01738 


0.17471 


0.26944 


At„ DD 


-0.09383 


-0.02507 


0.07411 


0.16636 


0.22319 




-0.09376 


-0.02506 


0.07405 


0.16624 


0.22304 


A c SD - no-cut 


-0.08977 


-0.02399 


0.07092 


0.15922 


0.21364 




-0.08968 


-0.02396 


0.07086 


0.15909 


0.21348 


A C FB SD - M 2 -cut 


-0.08978 


-0.02399 


0.07093 


0.15923 


0.21365 




-0.08968 


-0.02396 


0.07086 


0.15910 


0.21349 


A b DD 


0.03652 


0.06403 


0.10311 


0.13956 


0.16242 




0.03615 


0.06361 


0.10262 


0.13903 


0.16187 


A b FB SD - no-cut 


0.03556 


0.06233 


0.10035 


0.13579 


0.15802 




0.03529 


0.06209 


0.10014 


0.13562 


0.15788 


A b FB SD - M 2 -cut 


0.03556 


0.06233 


0.10035 


0.13580 


0.15803 




0.03529 


0.06209 


0.10014 


0.13562 


0.15788 



Table 14: T0PAZ0/ZFITTER comparison of A^ b in DD-mode and SD-modes. 



In Fig. we show the absolute deviation in per-mill between T0PAZ0 and 
ZFITTER predictions for the forward-backward asymmetries in DD mode and in 
SD-modes. From this figure one understands that A b B is indeed the RO showing 
the largest deviations between the two codes. This is hardly a surprise, given 
the comparison for sin 2 6 b s reported in Tab.(R). In both cases it is the absence 
of next-to-leading corrections in the 65-channel (due to missing non-universal 
next-to-leading terms) that stays at the root of the relatively large theoretical 
uncertainty. 

For A^ B the absolute deviations is always below 0.14 per-mill (0.04 — 0.05 
at the peak). For A FB the agreement is also very good, deviations below 0.16 
per-mill (and only at the wings) and peak asymmetries differing of 0.06 — 0.07 
per-mill. This sort of agreement and consistency between DD-mode and SD- 
modes shows that also final state QCD corrections are under control in the 
c-channel. 
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Figure 2: Absolute deviations between TDPAZO and ZFITTER for de-convoluted 
forward-backward asymmetries. 
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Note that QCD corrections for the forward-backward asymmetry can well 
be approximated by an expansion in the parameter x = Ami/ s = 4r 2 , 



•'va 



4a(i + ^/i). (53) 



with a correction factor, /i, which we write as [B6| 

8l/7 7T 2 2, 1, 2 \ 40 3 /55 7T 2 19, 1, 2 , ., 



3 ^3 18 3 3 V 27 \ 24 12 12 2 % 

(54) 
where l r — lnr. These first order corrections vanish in the massless limit. The 
asymmetry changes as 



A q ^ A q 



l-^(l-/i) . (55) 



For 6-quarks the inclusion of QCD final state correction improves the T0PAZ0- 
ZFITTER agreement. This fact is not completely satisfactory, signalling some 
difference (and some uncertainty) in the implementation of electroweak/QCD 
radiative corrections for the 66-channel. 

5.3 Higgs-Mass Dependence of De-Convoluted Observables 

Our comparison for de-convoluted ROs has to be extended to a wide range of 
values for the Higgs boson mass: differences indicate theoretical uncertainties in 
using the measurements to constrain the mass of the Higgs boson. In Tabs. ([15]- 
|l7|) we present cross-sections and asymmetries in the Higgs-mass range M H 
range of 10 — 1000 GeV. In Fig. g we show the relative deviations between the 
TDPAZO and ZFITTER prediction for a{fi), cr(had) and A^ B in SD no-cut mode 
as a function of the Higgs mass ranging from 10 GeV to 1000 GeV. 

The figure confirms the good agreement between the two sets of predictions: 
differences in the peak muonic cross-sections are everywhere below 0.14 per- 
mill, reached only at the boundaries of the interval in M H (0.05 per-mill at 
M H = 100 GeV). Deviations in the peak hadronic cross-sections vary from —0.13 
per-mill at very low values of the Higgs boson mass to —0.07 per-mill at M H = 
100 GeV and stay practically constant for higher values of M H . The variations 
in relative differences for <r(//) are 0.09 per-mill at peak and 0.16, 0.11 per mill at 
the wings (M z ± 3). For Ohad we have 0.07 per-mill at peak and 0.17, 0.10 per- 
mill at the wings. For A^ g we observe absolute differences which are everywhere 
below 0.00012 and at the peak below 0.00006. 
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a(li) [nb] 


V^[GeV] 


M H in GeV 


10 


30 


100 


300 


1000 


M z -'i 


0.30000 
0.30001 


0.30055 
0.30059 


0.30055 
0.30058 


0.30011 
0.30011 


0.29938 
0.29937 


M z - 1.8 


0.65724 
0.65726 


0.65832 
0.65839 


0.65839 
0.65845 


0.65766 
0.65766 


0.65643 
0.65641 


M z 


2.00558 
2.00537 


2.00613 
2.00601 


2.00711 
2.00700 


2.00827 
2.00809 


0.20098 
0.20095 


M z + 1.8 


0.65810 
0.65815 


0.65968 
0.65976 


0.65978 
0.65985 


0.65900 
0.65902 


0.65772 
0.65773 


M z +3 


0.31003 
0.31007 


0.31102 
0.31107 


0.31104 
0.31108 


0.31053 
0.31054 


0.30971 
0.30973 



Table 15: TOPAZO/ZFITTER comparison for the muonic cross-section in SD fully 
extrapolated mode for different values of the Higgs boson mass. 



cr(had) [nb] 


Vi[GeV] 


M^in GeV 


10 


30 


100 


300 


1000 


M z -3 


5.99467 
5.99740 


6.00515 
6.00774 


6.00291 
6.00518 


5.99098 
5.99267 


5.97237 
5.97414 


M z - 1.8 


13.40951 
13.41443 


13.42945 
13.43405 


13.42550 
13.42948 


13.40383 
13.40683 


13.37020 
13.37339 


M z 


41.42574 
41.43106 


41.42837 
41.43209 


41.43114 
41.43401 


41.43364 
41.43645 


41.43777 
41.44034 


M z + 1.8 


13.48804 
13.49074 


13.51762 
13.52007 


13.51353 
13.51559 


13.49007 
13.49146 


13.45456 
13.45640 


M z +3 


6.26327 
6.26441 


6.28209 
6.28316 


6.27961 
6.28050 


6.26553 
6.26606 


6.24431 
6.24513 



Table 16: TOPAZO/ZFITTER comparison for the hadronic cross-section in SD 
fully extrapolated mode for different values of the Higgs boson mass. 
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Figure 3: Relative deviations between TOPAZO and ZFITTER for cr{p), A^ B and 
er(had) for different values of the Higgs boson mass in SD no-cut mode. 
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A? 

FB 


^[GcV] 


M^in GeV 


10 


30 


100 


300 


1000 


M z -3 


-0.25984 


-0.26021 


-0.26122 


-0.26246 


-0.26398 




-0.25986 


-0.26019 


-0.26119 


-0.26246 


-0.26399 


M z - 1.8 


-0.14864 


-0.14909 


-0.15010 


-0.15126 


-0.15264 




-0.14868 


-0.14910 


-0.15010 


-0.15127 


-0.15265 


M z 


0.01899 


0.01838 


0.01742 


0.01644 


0.01534 




0.01893 


0.01832 


0.01737 


0.01641 


0.01533 


M z + 1.8 


0.17648 


0.17566 


0.17478 


0.17401 


0.17324 




0.17639 


0.17557 


0.17469 


0.17396 


0.17321 


M z +3 


0.27130 


0.27035 


0.26952 


0.26890 


0.26833 




0.27120 


0.27023 


0.26941 


0.26884 


0.26830 



Table 17: TOPAZO/ZFITTER comparison for the muonic forward-backward asym- 
metry in SD fully extrapolated mode for different values of the Higgs boson 
mass. 



5.4 Standard Model Remnants 

POs are determined by fitting ROs, but actually something is still taken from 
the SM (imaginary parts, parts which have been moved to interference terms 
and photon-exchange terms) making the model-independent results dependent 
on the SM. How complicated is such a description? Within the codes we consider 
sub-contributions to the DD de-convoluted quantities, 1) total DD, 2) DD with 
Z exchange only, 3) DD with Z © 7 without interference. For instance we 
construct the relative and absolute differences 



S int a 



r DD 



r DD 
,DD 



1 



in percent, 



,DD\ 



A WB [ A FB ) z-\ 



(56) 



They are reported in Tab.(|18[). The effect of the Z — 7 interference is negative 
below the peak, vanishingly small around it and turning positive and large above 
it. The effect of Z — 7 interference is particularly important for the forward- 
backward asymmetry, as its energy dependence is governed by this interference. 

Among the de-convoluted quantities the most relevant are those computed 
at s = M 2 , which have an obvious counterpart in the PO, that we have already 
computed, i.e., a\ , a\, and A°£. 

There is however a noticeable difference between the two sets, represented 
by the interference of the Z — 7 s-channel diagrams, including the imaginary 
parts in a(s) and in the form-factors, the latter being particularly relevant for 
the leptonic forward-backward asymmetry. This effect is illustrated in Tab. (19). 
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Centre-of-mass energy in GeV 


M z - 3 


M z - 1.8 


M z 


M z + 1.8 


M z + 3 


6 [nt a^ 


-0.209% 


-0.136% 


-0.028% 


+0.072 % 


+0.132% 


<5 in V had 


-0.492 % 


-0.301 % 


-0.029% 


+0.226 % 


+0.384% 


A int A^ 

FB 


-0.27703 


-0.16611 


0.00145 


0.15923 


0.25441 



Table 18: T0PAZ0 relative differences S lnt a and absolute differences A lnt A£ B . 





Centre-of-mass energy in GeV 


M z - 3 


M z - 1.8 


M z 


M z + 1.8 


M z + 3 


(T M [nb] No Ims 
cr M [nb] 


0.29996 
0.29999 


0.65713 
0.65718 


2.00343 
2.00341 


0.65855 
0.65856 


0.31045 
0.31047 


Diff.[pb] 


+0.03 


+0.05 


-0.02 


-0.01 


+0.02 


Chad [nb] No Ims 
Chad [nb] 


5.78583 
5.78492 


12.94061 
12.93841 


39.93848 
39.92967 


13.02635 
13.02421 


6.05322 
6.05233 


Diff.[pb] 


-0.91 


-2.20 


-8.81 


-2.14 


-0.89 


A^ No Ims 

FB 

A>" 

FB 


-0.26311 
-0.26170 


-0.15181 
-0.15037 


0.01598 
0.01745 


0.17364 
0.17510 


0.26858 
0.27002 


Diff. 


-0.00141 


-0.00144 


-0.00147 


+0.00146 


-0.00144 



Table 19: T0PAZ0 comparison of DD (completely de-convoluted) RO 

with/without imaginary parts (Ims) in couplings and form factors. 



5.5 The Z — 7 Interference for Cross-Sections 

One must evaluate the residual SM dependence of the so-called model-inde- 
pendent parameters; one simple source of such a SM dependence is due to the 
interference terms for cross-sections, which are governed by the value of sin 6*g 
and therefore depend on the values of nit and M H chosen.g Note that for leptonic 
final states, one can use the POs to express the interference terms, at least up 
to imaginary parts which must be taken from the SM as usual. This is possible 
because the interference terms are proportional to the effective couplings which 
can be derived from the POs Re and -Apg-n However, for the inclusive hadronic 
final state, which is a sum over the five light quark flavours, the interference 
terms must be taken completely from the SM.n 



2 The LEP community has agreed on a set of numbers, mt = 175 GeV and M H = 150 GeV. 

3 This is realised for MI calculations with TOPAZO. For MI calculations with ZFITTER, it is 
realised for the effective-couplings interfaces ,_hutnot for the partial-width interface. 

4 This is avoided in the S-Matrix ansatz [Eq, b9| , which treats also the Z — 7 interference 
terms for cross sections and asymmetries as free and independent parameters to be determined 
from the data. The experimental measurements are also analysed within this extended MI 
ansatz. Combined LEP results are given in ||13|. 
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M H [GcV] 


Centre-of-mass energy in GeV 


M z - 3 


M z - 1.8 


M z 


M z + 1.8 


M z + 3 


S'^a^ 


10 
100 
1000 


-0.229% 
-0.209% 
-0.181% 


-0.150% 
-0.136% 
-0.119% 


-0.030% 
-0.028% 
-0.026% 


+0.082 % 
+0.072 % 
+0.060 % 


+0.149% 
+0.132% 
+0.111% 


<s in v had 


10 
100 
1000 


-0.518% 
-0.492 % 
-0.457 % 


-0.317% 
-0.301 % 
-0.281 % 


-0.029% 
-0.029% 
-0.028% 


+0.240 % 
+0.226 % 
+0.207% 


+0.407% 
+0.384% 
+0.353% 



Table 20: T0PAZ0 relative differences S int a(M H ) = o- DD /cr^ 7 - 1 as a function 
of the Higgs boson mass in DD-mode. 



In Tab. (20) we show the relative deviation of excluding/including the Z — 7 
interference as a function of the Higgs boson mass in DD-mode. As observed 
before the Z — 7 interference is negative below the peak and changes sign above 
it. It is vanishingly small at the resonance for all values of M H , approximately 
—0.03%, and can be sizeable at the wings, — 0.2%(— 0.5%) at the left wing and 
+0.1%(+0.4%) at the right wing for the muonic (hadronic) cross-section. The 
Higgs-mass dependence of the Z — 7 interference is rather large, up to 20% of 
the interference itself. 



6 Convoluted Realistic Observables 

6.1 Comparison for Extrapolated Setup 

Having discussed the status of our comparisons before the introduction of initial- 
state QED radiation we now proceed to comparing the convoluted quantities 
and the effect of convolution. 

The default of T0PAZ0/ZFITTER is to account for initial-state QED radiation 
through a so-called additive formulation of the QED radiator ( flux- function) , 
which is a mixture of leading-logarithms (LL) and finite-order results. In J3(J 
a proof is given that the /3a;' 3-1 term should be factorized in front of, at least, 
the LL component. This result is obtained up to third order LL, and there are 
good indications from fourth and fifth orders that it is true to infinite order.^ 
Recently explicit a 3 L 3 terms became known |}0) and also |31j] . For higher orders 
we refer to [ J32J an d J33|. In TOPAZO/ZFITTERthe O (a 3 ) radiator is implemented 
according to Q. Recently TDPAZO and ZFITTER have implemented the order 
a 3 factorized (YFS) radiator as reported in |3J| . 

5 S. Jadach, private communication. 
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There is a pattern of convolution that we want to compare in our step-by- 
step procedure. 

CA3 Complete RO, with QED initial state radiation implemented through an 
additive O (a 3 ) radiator Q. 

CF3 Complete RO, with QED initial state radiation implemented through a 
factorized O (a 3 ) radiator fj5| . 

The following two equations define cross-sections and forward-backward asym- 
metries convoluted with ISR: 



/ dzH (z;s) a T (zs), 



a T (s) = / dzH (z; s) <t t (zs), (57) 



where z = sq/s and 



ira 2 Q 2 e Q 2 f [ l 1 . , 

4b (s) = / dz- -g H FB (z; s) cr pB (zs) , (58) 

°tot Jz (1 + Z) 

Note that the so-called radiator (or flux function) , H (z; s), is known up to terms 
of order a 3 while H FB is only known up to terms of order a 2 . 

The kernel cross-sections <r T FB should be understood as the improved Born 
approximation (IB A), including imaginary parts and corrected with all elec- 
troweak and possibly all FSR (QED (g) QCD) corrections where all coupling 
constants and effective vector and axial weak couplings are running, i.e., they 
depend on s' — zs under the convolution integrals in Eq.(j57|) and Eq.(pq). In 
practice this takes a lot of CPU time and for this reason some time-saving op- 
tions are foreseen in the codes. For instance, one may calculate effective weak 
couplings only once at s rather than at s' thereby saving a conspicuous amount 
of CPU time. 

We study the accuracy of such approximations with ZFITTER. In Tab. |T] 
we report ROs calculated with no convolution at all (all couplings evaluated 
at s), convolution of a only (a(s) — > a(s')), and full convolution of all elec- 
troweak radiative corrections; corresponding to the ZFITTER flag CONV with 
values C0NV=0,1,2, respectively. The bulk of the running-couplings effect is 
given by the a convolution, in particular below the wing. The remaining effect 
is totally negligible at the resonance and below, growing to ~ —0.15 per- mill 
above the resonance for cross-sections and remaining negligible for the asymme- 
tries. This study proves that one may avoid using the CPU-time consuming full 
convolution of all electroweak radiative corrections and that it is sufficient to 
keep the a convolution only. This is the default used for the ZFITTER numbers 
reported in this article. 

In T0PAZ0 all universal electroweak corrections and final-state QCD correc- 
tions are put in convolution with initial-state QED radiation and therefore the 
couplings, a, and a s are evaluated at the scale s'. Weak boxes, vertices and 
expanded bosonic self-energy corrections are added linearly, evaluated at the 
nominal energy. The latter is also true for IFI. 
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Centre-of- 


mass energy in GeV 




M z -3 


M z - 1.8 


M z 


M z + 1.8 


M z +3 


o> [nb] 


0.22862 


0.47672 


1.48012 


0.69526 


0.40655 




0.22843 


0.47653 


1.47995 


0.69509 


0.40638 




0.22842 


0.47653 


1.47992 


0.69505 


0.40633 




-0.04 





-0.02 


-0.06 


-0.12 


Chad [nb] 


4.45219 


9.60235 


30.43892 


14.18454 


8.19986 




4.45146 


9.60165 


30.43824 


14.18391 


8.19923 




4.45130 


9.60139 


30.43753 


14.18270 


8.19791 




-0.04 


-0.03 


-0.02 


-0.09 


-0.16 


AL 


-0.28308 


-0.16979 


-0.00066 


0.11177 


0.15451 




-0.28330 


-0.16985 


-0.00066 


0.11182 


0.15461 




-0.28330 


-0.16985 


-0.00066 


0.11183 


0.15464 













+0.01 


+0.03 



Table 21: ZFITTER illustration of the effect of convolution of electroweak ra- 
diative corrections in CA3 mode with a cut of s' > 0.01s. First line is no 
convolution, second line is a(s) — > a(s') convolution, third line is full convolu- 
tion of all electroweak radiative corrections. Fourth line is difference between 
third and second line. 



Results are shown in Tab.(|22|) where we include all steps in the process de- 
convoluted — ► convoluted, i.e., DD, SD, CA3 and CF3. The reported results 
refer to an M 2 (//) cut of 0.01 s. The relative deviations between T0PAZ0 and 
ZFITTER are shown in Fig. [§ The numbers of Tab.@) are produced with a 
M 2 (//)-cut of 0.01 s. It is instructive to compare with similar results obtained 
by imposing the an s'-cut, e.g. s' > 0.01 s. The comparison is shown in Tab.(f23|) 
and Fig. ||. From Tab.(|23j) one sees that the differences between the two cuts are 
of 0.43 (1.38) pb for the muonic (hadronic) peak cross-section. There is really 
no problem as long as the procedure is fully specified. 

It emerges from this comparisons that the agreement for the muonic cross- 
section for energies below the peak is quite reasonable but less satisfactory than 
for hadrons. Given the agreement at the level of de-convoluted cross-sections 
and once we have observed that the de-convolution is satisfactory for hadrons, 
we come to the conclusion that for muons the low-q 2 region, where the Coulomb 
pole and mass effects (for very loose s'-cuts) may become relevant, gives the 
dominant difference in a a . 
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Centre-of- 


mass energy in GeV 




M z -3 


M z - 1.8 


M z 


M z + 1.8 


M z +3 


<V [ nb ] DD 


0.29999 


0.65718 


2.00341 


0.65856 


0.31047 




0.30003 


0.65724 


2.00331 


0.65863 


0.31051 


*>, [nb] SD 


0.30047 


0.65821 


2.00656 


0.65960 


0.31095 




0.30052 


0.65832 


2.00659 


0.65971 


0.31102 


cr M [nb] CA3 


0.22840 


0.47642 


1.47967 


0.69490 


0.40628 




0.22836 


0.47641 


1.47962 


0.69492 


0.40627 


c^ [nb] CF3 


0.22841 


0.47645 


1.47977 


0.69495 


0.40630 




0.22837 


0.47644 


1.47971 


0.69497 


0.40629 


Chad [nb] DD 


5.78492 


12.93841 


39.92967 


13.02421 


6.05233 




5.78670 


12.94148 


39.93079 


13.02591 


6.05313 


Chad [nb] SD 


6.00265 


13.42490 


41.42929 


13.51293 


6.27933 




6.00500 


13.42901 


41.43272 


13.51517 


6.28030 


o-had [nb] CA3 


4.44990 


9.59865 


30.43501 


14.18203 


8.19853 




4.45129 


9.60132 


30.43725 


14.18342 


8.19894 


(Thad [nb] CF3 


4.45016 


9.59921 


30.43696 


14.18307 


8.19901 




4.45157 


9.60191 


30.43929 


14.18451 


8.19945 


A» DD 


-0.26170 


-0.15037 


0.01745 


0.17510 


0.27002 




-0.26167 


-0.15037 


0.01741 


0.17502 


0.26991 


^SD 


-0.26128 


-0.15013 


0.01742 


0.17481 


0.26958 




-0.26122 


-0.15011 


0.01738 


0.17471 


0.26944 


A» B CAS 


-0.28321 


-0.16981 


-0.00062 


0.11189 


0.15470 




-0.28336 


-0.16988 


-0.00066 


0.11184 


0.15464 


,4£ B CF3 


-0.28320 


-0.16980 


-0.00062 


0.11189 


0.15469 




-0.28340 


-0.16990 


-0.00066 


0.11185 


0.15465 



Table 22: TOPAZO/ZFITTER comparison of complete RO (CA3 and CF3) with 
DD and SD (M 2 -cut) modes. 
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Figure 4: Relative deviations between TDPAZO and ZFITTER for muonic and 
hadronic cross-section in DD, SD, CA3 and CF3 modes. The last three are 
corresponding to a cut M 2 (ff) of 0.01 s. 
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Centre-of- 


mass energy in GeV 




M z -3 


M z - 1.8 


M z 


M z + 1.8 


M z +3 


<V [ nb ] DD 


0.29999 


0.65718 


2.00341 


0.65856 


0.31047 




0.30003 


0.65724 


2.00331 


0.65863 


0.31051 


*>, [nb] SD 


0.30055 


0.65839 


2.00711 


0.65978 


0.31104 




0.30058 


0.65844 


2.00700 


0.65985 


0.31108 


cr M [nb] CA3 


0.22849 


0.47657 


1.48010 


0.69512 


0.40642 




0.22843 


0.47653 


1.47995 


0.69509 


0.40638 


c^ [nb] CF3 


0.22850 


0.47660 


1.48019 


0.69517 


0.40644 




0.22844 


0.47656 


1.48004 


0.69514 


0.40640 


Chad [nb] DD 


5.78492 


12.93841 


39.92967 


13.02421 


6.05233 




5.78670 


12.94148 


39.93079 


13.02591 


6.05313 


Chad [nb] SD 


6.00291 


13.42550 


41.43114 


13.51353 


6.27961 




6.00518 


13.42948 


41.43401 


13.51559 


6.28050 


o-had [nb] CA3 


4.45012 


9.59910 


30.43639 


14.18269 


8.19892 




4.45146 


9.60165 


30.43824 


14.18391 


8.19923 


(Thad [nb] CF3 


4.45038 


9.59966 


30.43834 


14.18373 


8.19940 




4.45174 


9.60225 


30.44028 


14.18499 


8.19974 


A» DD 


-0.26170 


-0.15037 


0.01745 


0.17510 


0.27002 




-0.26167 


-0.15037 


0.01741 


0.17502 


0.26991 


^SD 


-0.26122 


-0.15010 


0.01742 


0.17478 


0.26952 




-0.26119 


-0.15010 


0.01737 


0.17469 


0.26941 


A» B CAS 


-0.28312 


-0.16977 


-0.00062 


0.11186 


0.15466 




-0.28330 


-0.16985 


-0.00066 


0.11182 


0.15461 


,4£ B CF3 


-0.28311 


-0.16976 


-0.00062 


0.11186 


0.15465 




-0.28333 


-0.16987 


-0.00066 


0.11183 


0.15462 



Table 23: TOPAZO/ZFITTER comparison of complete RO (CA3 and CF3) with 
DD and SD (s'-cut) modes. 
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Figure 5: Relative deviations between TDPAZO and ZFITTER for muonic and 
hadronic cross-section in DD, SD, CA3 and CF3 modes. The last three are 
corresponding to a cut on s' of 0.01 s. 
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The effect of moving the s'-cut from 0.01 s to 0.1 s is shown in Tab. (124). As 
a result the agreement becomes much better, especially on the lower tail of the 
resonance. This fact is also shown in Fig. pi 







Centre-of- 


mass energy in GeV 




M z - 3 


M z - 1.8 


M z 


M z + 1.8 


M z + 3 


* M [lib] SD 


0.30055 


0.65839 


2.00711 


0.65978 


0.31104 


cr M [nb] CA3 


0.30058 


0.65845 


2.00700 


0.65985 


0.31108 


0.22849 


0.47657 


1.48010 


0.69512 


0.40642 




0.22843 


0.47653 


1.47995 


0.69509 


0.40638 




0.22674 


0.47487 


1.47845 


0.69353 


0.40487 




0.22674 


0.47489 


1.47836 


0.69356 


0.40489 


^sd 


-0.26122 


-0.15010 


0.01742 


0.17478 


0.26952 


A^ CAS 


-0.26119 


-0.15010 


0.01737 


0.17469 


0.26941 


-0.28312 


-0.16977 


-0.00062 


0.11186 


0.15466 




-0.28330 


-0.16985 


-0.00066 


0.11182 


0.15461 




-0.28526 


-0.17035 


-0.00061 


0.11214 


0.15528 




-0.28532 


-0.17040 


-0.00064 


0.11210 


0.15523 



Table 24: T0PAZ0/ZFITTER comparison of complete RO (CA3-mode) with SD 
(s'-cut) mode. For CA3-mode first (second) row correspond to T(Z) s' > 0.01 s. 
Third (fourth) row give instead s' > 0.1 s. 



In Tab. (25) we show the comparison for heavy quark forward-backward 



asymmetries, including initial-state QED radiation. As we have observed before 
the ^-channel shows larger deviations. Even though the agreement for the con- 
voluted A b is satisfactory, especially at the peak, one should not forget that a 
similar comparison for DD-de- convoluted A^ B is considerably worse so that the 
result of Tab. (25) is also a consequence of accidental compensations. 





Centre-of-mass energy in GeV 


M z -3 


M z - 1.8 


M z 


M z + 1.8 


M z +3 


K» 


-0.10600 

-0.10598 

-0.02 


-0.03625 

-0.03625 

0.00 


0.06068 

0.06065 

+0.03 


0.12386 

0.12377 

+0.09 


0.14840 

0.14827 

+0.13 


< 


0.028131 

0.028078 

+0.05 


0.05705 

0.05701 

+0.04 


0.09611 

0.09612 

-0.01 


0.12135 

0.12161 

-0.26 


0.13105 

0.13169 

-0.64 



Table 25: T0PAZ0/ZFITTER comparison of heavy quark forward-backward 
asymmetries in CA3-mode. Third row is the absolute deviation in per-mill. 



47 




0.2 



— . 0.1 



£ 0.0 - 
Q. 



N 

I- 



■0.1 



"l — n r 



T~ 



SD s* = 0.01 s 



SD s' - 0.1 s 



DD 



-0.2 

87 88 89 90 91 92 93 94 95 

E rm (GeV) 

cm 

Figure 6: Relative deviations between TDPAZO and ZFITTER for muonic to- 
tal cross-section and forward-backward asymmetry in SD-mode and CA3-mode 
with two set of s'-cuts. 
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6.2 Comparison with Realistic Kinematical Cuts 

We have also devoted an effort in order to present the most up-to-date analysis 
for ROs with realistic kinematical cuts. Therefore, we go beyond the fully 
extrapolated set-up for muonic channel (with the inclusion of a s' /M 2 -cut) by 
considering 



H+H' 



for0 a 



< 0- < 7T • 



0°,20° and 40°) 



3 il < 10°, 25°, and EthifJ^) = 1 GeV, 



where 9- is the final-state fermion scattering angle and # aco n the acollinearity 
between the final-state fermions. The results are shown in Tabs.(126H2T 





^acc 




Ccntre-of- 


mass energy in GeV 




M z - 3 


M z - 1.8 


M z 


M z + 1.8 


M z + 3 


o> [nb] 


0° 


0.21932 


0.46287 


1.44795 


0.67725 


0.39366 






0.21928 


0.46285 


1.44781 


0.67722 


0.39361 






+0.18 


+0.04 


+0.10 


+0.04 


+0.13 


20° 


0.19990 


0.42207 


1.32066 


0.61759 


0.35886 






0.19987 


0.42205 


1.32053 


0.61756 


0.35881 






+0.15 


+0.05 


+0.10 


+0.05 


+0.14 


40° 


0.15034 


0.31762 


0.99428 


0.46479 


0.26989 






0.15032 


0.31760 


0.99415 


0.46474 


0.26983 






+0.13 


+0.06 


+0.13 


+0.11 


+0.22 


AL 


0° 


-0.28450 


-0.16914 


0.00033 


0.11512 


0.16107 






-0.28453 


-0.16911 


0.00025 


0.11486 


0.16071 






+0.03 


-0.03 


+0.08 


+0.26 


+0.36 


20° 


-0.27509 


-0.16352 


0.00042 


0.11171 


0.15645 






-0.27506 


-0.16347 


0.00035 


0.11148 


0.15616 






-0.03 


-0.05 


+0.07 


+0.23 


+0.29 


40° 


-0.24219 


-0.14396 


0.00054 


0.09906 


0.13903 






-0.24207 


-0.14386 


0.00050 


0.09893 


0.13891 






-0.12 


-0.10 


+0.04 


+0.13 


+0.12 



Table 26: T0PAZ0/ZFITTER comparison for muonic total cross-section and 
forward-backward asymmetry of complete RO (CA3-mode) with the angular 
acceptance (9 acc = 0°,10°,20°) and acollinearity (9 aco \ < 10°) cuts. First row 
T0PAZ0, second row ZFITTER, third row relative (absolute) deviations in per- 
mill. 

We register an agreement comparable with the one obtained with s' /M 2 - 
cut, perhaps deteriorating a little for A^ at the wings. In conclusion the 
agreement between T0PAZ0 and ZFITTER remains rather remarkable even when 
the geometrical acceptance is constrained and also final-state energies and the 
acollinearity angle are bounded. 
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"ace 




Centre-of- 


mass energy in GeV 




M z -3 


M z - 1.8 


M z 


M z + 1.8 


M z +3 


o> [iib] 


0° 


0.22333 


0.46971 


1.46611 


0.68690 


0.40034 






0.22328 


0.46968 


1.46598 


0.68688 


0.40031 






+0.22 


+0.06 


+0.09 


+0.03 


+0.075 


20° 


0.20359 


0.42835 


1.33731 


0.62648 


0.36507 






0.20357 


0.42833 


1.33718 


0.62647 


0.36505 






+0.10 


+0.05 


+0.10 


+0.02 


+0.055 


40° 


0.15320 


0.32245 


1.00698 


0.47167 


0.27479 






0.15318 


0.32243 


1.00682 


0.47164 


0.27477 






+0.13 


+0.06 


+0.16 


+0.06 


+0.07 


A£ 


0° 


-0.28617 


-0.17037 


-0.00032 


0.11324 


0.15730 






-0.28647 


-0.17049 


-0.00043 


0.11293 


0.15682 






+0.30 


+0.12 


+0.11 


+0.31 


+0.48 


20° 


-0.27695 


-0.16485 


-0.00026 


0.10974 


0.15250 






-0.27722 


-0.16497 


-0.00037 


0.10944 


0.15204 






+0.27 


+0.12 


+0.11 


+0.30 


+0.46 


40° 


-0.24423 


-0.14536 


-0.00016 


0.09703 


0.13492 






-0.24445 


-0.14545 


-0.00026 


0.09678 


0.13454 






+0.22 


+0.09 


+0.10 


+0.25 


+0.38 



Table 27: The same as in Tab. 26 but for the acollinearity cut 6> aco i < 25° 



We note that the coding in ZFITTER, for the part involving realistic cuts, is 
based on some old work [[56J . A recent study, presented in Q , shows that the 
approximations made in the former reference ensure sufficient technical precision 
of the treatment of ISR, O (l0~ 4 ), at SLD/LEP-1 energies. (See Section g for 
the situation concerning initial-final QED interference). Coding in TDPAZO is 
always based on the work of p3{ . 

6.3 Uncertainty on QED Convolution 

We now return to a detailed analysis of initial-state QED radiation by defining 
convolution factors for each realistic observable O, giving the net effect of initial- 
state QED radiation at the various energies. 



edec 



(O) 



o 

QSD 



1, 



A dec (0) = 0-0 



SD 



(59) 



For convenience of the reader we reproduce in Tab.(|28|) the results for the CA3 
and CF3 mode. From Tab.(Pq) we derive the absolute differences, for A 11 , and 
the relative ones, for cross-sections, between additive and factorized versions of 
the QED radiators. They are shown in Tab. (129). 
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Centre-of-mass energy in GeV 


M z -3 


M z - 1.8 


M z 


M z + 1.8 


M z +3 


£dec( a CA3) T 

S dcc (a^ A3 ) Z 
jdcc^CFS) T 

jdcc^CFS) z 


-23.976 
-24.007 
-23.973 
-24.000 


-27.616 
-27.629 
-27.611 
-27.624 


-26.257 
-26.261 
-26.253 
-26.256 


5.356 
5.341 
5.364 
5.348 


30.665 
30.631 
30.671 
30.637 


^ 0C K C a A d 3 ) T 
^c K CA 3) z 

^ CC (^a F d 3 ) Z 


-25.867 
-25.873 
-25.863 
-25.869 


-28.501 
-28.503 
-28.497 
-28.499 


-26.537 
-26.538 
-26.533 
-26.533 


4.952 
4.945 
4.959 
4.953 


30.564 
30.550 
30.572 
30.558 


^dcc/^CA3-\ T 
A dcc (A^A3) z 
^dec/^CF3\ T 
A dcc (A^F3) z 


-2.190 
-2.211 
-2.189 
-2.215 


-1.967 
-1.975 
-1.966 
-1.977 


-1.804 
-1.803 
-1.804 
-1.803 


-6.292 
-6.287 
-6.292 
-6.286 


-11.486 
-11.480 
-11.487 
-11.479 



Table 28: 
modes. 



RO: the effect in % of initial state QED radiation for CA3 and CF3 





Centre-of-mass energy in GeV 


M z -3 


M z -\ 


M z 


M z + 1 


M z +3 


10 4 x (fact/add-1) 


cv 


0.44 
0.88 


0.63 
0.63 


0.61 
0.68 


0.72 
0.72 


0.49 
0.49 


Chad 


0.58 
0.61 


0.58 
0.62 


0.64 
0.67 


0.73 
0.76 


0.59 
0.62 


fact-add [pb] 


o> 


0.01 
0.02 


0.03 
0.03 


0.09 
0.10 


0.05 
0.05 


0.02 
0.02 


Chad 


0.26 
0.27 


0.56 
0.60 


1.95 
2.04 


1.04 
1.08 


0.48 
0.51 


10 5 x (fact-add) 


A* 

FB 


1.00 
-4.00 


1.00 
-2.00 


0.00 
0.00 


0.00 
1.00 


-1.00 
1.00 



Table 29: Absolute and relative differences in TDPAZO and in ZFITTER for 
additive and factorized radiators. 
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To give an example of the developments in the treatment of QED initial-state 
radiation we recall that in TDPAZO the following steps have occurred: 

• the leading O (a 2 ) result was considered in version 1.0, 

• the complete O (a 2 ) result, i.e., leading plus NLO O (a 2 L) and NNLO 
O (a 2 L°) was added in version 2.0 ||§, 

• the complete O (a 2 ) plus leading O (a 3 !/ 3 ) result as been included after 
version 4.0 ||, 

where L = ln(s/m 2 ). Inserting the O (a 3 L 3 ) terms into the additive radiator 
leads to a negative shift of —0.59 per-mill in the peak hadronic cross-section. 
The complete shift leading-a 2 to leading-a 3 is dominated by the leading O (a 3 ) 
terms with very little influence by the NLO O (a 2 L) terms. 

Sometimes the forward and backward cross-sections rx F B , are actually used 
to calculate ^ FB - In Tab.(|30|) we present results for the total/F/B cross-sections 
obtained with the additive radiator (CA3), showing the effect due to initial-state 
radiation for the forward and backward cross-section separately. Comparing 
results obtained with the additive and factorized radiators, we estimate the 
corresponding initial-state QED uncertainty as reported in Tab. all. 



a [nb] 




Centre-of- 


mass ener 


gy in GeV 




M z -3 


M z - 1.8 


M z 


M z + 1.8 


M z +3 


CT M 


0.22849 


0.47657 


1.48010 


0.69512 


0.40642 


< 


0.08190 


0.19783 


0.73959 


0.38644 


0.23464 




-26.229 


-26.661 


-26.707 


-22.655 


-25.570 


< 


0.14659 


0.27874 


0.74050 


0.30868 


0.17178 




-22.655 


-25.570 


-24.260 


+13.388 


+51.210 


Gb 


0.95245 


2.06397 


6.56297 


3.05537 


1.76374 


< 


0.48962 


1.09086 


3.59687 


1.71307 


0.99744 




-29.290 


-28.973 


-29.175 


-26.377 


-28.229 


*? 


0.46283 


0.97311 


2.96610 


1.34229 


0.76630 




-26.377 


-28.229 


-27.413 


+6.642 


+34.719 


<Jc 


0.77811 


1.66456 


5.25261 


2.45308 


1.42235 


^ 


0.34781 


0.80211 


2.78568 


1.37846 


0.81671 




-27.564 


-26.845 


-27.176 


-24.904 


-26.237 


«? 


0.43029 


0.86245 


2.46693 


1.07461 


0.60563 




-24.904 


-26.237 


-25.666 


+9.428 


+41.341 



Table 30: Total, forward and backward a f from T0PAZ0 (CA3 mode): first 
entry is the cross-section, second entry is the effect in % of the initial state 
QED convolution. 
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<x[pb] 


Centre-of-mass energy in GeV 


M z - 3 


M z - 1.8 


M z 


M z + 1.8 


M z +3 


°n 


0.01 


0.03 


0.09 


0.05 


0.02 




0.02 


0.03 


0.10 


0.05 


0.02 


°l 


0.01 


0.01 


0.05 


0.02 


0.01 




0.00 


0.01 


0.05 


0.03 


0.02 


< 


0.01 


0.01 


0.05 


0.02 


0.01 




0.01 


0.02 


0.05 


0.02 


0.01 


°b 


0.06 


0.12 


0.42 


0.22 


0.10 




0.06 


0.13 


0.44 


0.24 


0.11 


<tf 


0.03 


0.06 


0.21 


0.11 


0.05 




0.03 


0.07 


0.28 


0.15 


0.07 


*f 


0.03 


0.06 


0.21 


0.11 


0.05 




0.03 


0.05 


0.17 


0.09 


0.04 


Oc 


0.04 


0.10 


0.33 


0.18 


0.08 




0.05 


0.10 


0.35 


0.19 


0.09 


«£ 


0.02 


0.05 


0.17 


0.09 


0.04 




0.01 


0.05 


0.21 


0.13 


0.06 


*? 


0.02 


0.05 


0.17 


0.09 


0.04 




0.03 


0.06 


0.15 


0.06 


0.02 



Table 31: QED IS uncertainties, CF3-CA3, in pb, for total, forward and 
backward 07 from T0PAZ0, ZFITTER. 



6.4 Higgs-Mass Dependence of Convoluted Observables 

Comparisons for Higgs masses in the range from 10 GeV to 1 TeV are needed. 
For the error determination on M H in SM fits, for example, the M H varia- 
tion extends to that region, thus the calculation must be reliable there. Also, 
in principle, all the other SM parameters should be varied to make sure the 
(dis)agreement is not too dependent on actual central values. In Fig. |^ we 
repeat a comparison for different values of the Higgs boson mass, this time 
including initial-state QED radiation in CA3 mode with an s'-cut of 0.01 s. 

For the hadronic cross-section the relative difference T0PAZ0 - ZFITTER varies 
of 0.07 per-mill at the peak and over the whole interval 10 GeV — 1 TeV being 
-0.06 per-mill at M H — 100 GeV. At the boundaries of the interval in M H the 
difference register a variation of 0.18 per-mill at M H — 10 GeV (being -0.29 
per-mill at M H = 100 GeV) and of 0.09 per mill at M H = 1 TeV (being -0.04 
per mill at M H = 100 GeV). 

The results show that there is a tiny M H -dependence in the T0PAZ0-ZFITTER 
comparison at a fixed energy and that the differences are comparatively larger 
below the peak and for very low or very large values of M H . 
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Figure 7: Relative deviations between TDPAZO and ZFITTER for hadronic cross- 
section in CA3 mode as a function of the Higgs boson mass. The curves corre- 
spond to a cut s' > 0.01 s. 
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7 Initial-Final QED Interference 

If initial-final QED (ISR-FSR) interference (IFI) is included in a calculation 
there is a conceptual problem with the meaning of the s'-cut. In this case the 
definition of the variable s' is unnatural since one does no longer know the origin 
of the radiative photon (ISR or FSR). Only a cut on the invariant mass of the 
final-state // system, M(ff), makes sense. 

There is another option: to select events with little initial-state radiation 
one can use a cut on the acollinearity angle # aco i, between the outgoing fcrmion 
and anti-fermion. A cut on # aco i is roughly equivalent to a cut on the invariant 
mass of the // system, indeed one may write 

Af_ 1 - sin(fl acol /2) 

a ~ Zcff_ l + sin(0 acol /2)' [W) 

therefore a cut of # aco i < 10° is roughly corresponding to the request that 
M 2 /s > 0.84. 

The inclusion of initial-final QED interference in T0PAZ0 and in ZFITTER 
is done at O (a), i.e., the O (a) initial-final interference term is added linearly 
to the cross-sections without entering the convolution with ISR and without 
cross-talk to FSR. For loose cuts the induced uncertainty is rather small and 
the effect of the interference itself is minute, as shown below. On the contrary, 
when we select a tight acollinearity cut the resulting limit on the energy of the 
emitted photon becomes more stringent and the effect of interference grows. The 
corresponding theoretical uncertainty, due to missing higher-order corrections, 
is therefore expected to be larger. 

A pragmatic point of view to improve upon the pure O (a) inclusion of 
the interference, which however is not implemented in the codes, would be the 
following: let us consider corrections not belonging to the two classes of ISR and 
FSR. These non-factorizable corrections correspond to both interference terms 
and QED box diagrams. We have, therefore, three classes of contributions for 
the total cross-section: ISR, FSR and interference. 

If we denote by a K the sum of factorizable (Born) plus non-factorizable cross- 
sections the totally radiatively corrected process can be effectively described in 
terms of structure functions. In this case s' is still the square of the invari- 
ant mass available after initial-state QED radiation and s' > M 2 (ff) > So. 
Therefore we can write 

a{s)= I dzH in {z, S )H fin {z-z 01 s)a K (s;M 2 (f7)> So ), (61) 

Jz o 

where i?i n ,fin is the initial- and final-state QED radiator and z — M 2 (ff)/s. 
Moreover a K is the kernel cross-section evaluated with the constraint M 2 (ff) > 
So and s — zs. The accuracy obtainable with Eq.(pTI) depends on the various 
contributions included within the three factors, the initial-state radiator H m , 
the final-state radiator H^ n and the kernel cross-section including factorizable 
and non-factorizable parts. 
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7.1 Comparison for Extrapolated Setup 

The effect of including or excluding initial-final state QED interference is shown 
in Tab.(p2[). As seen from the Table, the level of agreement between TDPAZO 
and ZFITTER does not deteriorate after the inclusion of initial-final state QED 
interference for loose M 2 -cuts. For this setup the effect of O (a) IFI is under 
control. 6 

We also perform other comparisons to understand the effect of IFI for various 
setups. First of all we use an extrapolated setup with a tight M 2 -cut of 0.8 s 
(it was 0.01 s before). Note that such a cut is, not exactly, but approximately 
equivalent to # aco i = 10° (note that an M 2 -cut takes into account only the 
isotropic part of the photon phase space). The result, reported in Tab.(|33|), 
shows that the two codes nicely agree for the net IFI effect also in case of tight 
M 2 -cuts. Note that the effect of QED initial-final interference is quite sizeable, 
differently from what happens with loose cuts. 

A further test, performed with TOPAZO, selects another setup. Note that 
E th = 1 GeV at ^fs = M z - 3 is equivalent to an M 2 -cut of M 2 (ff) « 
2 E t h\/~s ~ 13.83997 GeV if we consider the radiative photon collinear with one 
of the final-state fermions. In this example we apply an M 2 -cut as above and 
a cut of 10° on the acollinearity so that the photon phase space is not isotropic 
anymore. The result is shown in Tab.(p4[). 

The cut on acollinearity, superimposed on the M 2 -cut, hardly changes the 

size of IFI effect, at least as seen by TOPAZO. Roughly speaking, one computes 

IFI with a cut on M 2 simulated by the acollinearity and has to add a small 

portion of the photon phase space, subtracting, at the same time, the region 

not allowed by acol-cut. To understand which portion of photon phase-space is 

added we introduce variables 

1 , . _ 1 _ 1 

- (1 - xi + x 2 ) Vs , E- = -xiVs , E i=2 

and observe that this parametrisation has the advantage that in the limit where 
we neglect mj the boundaries of the phase space are of a triangular form 

%i = x-i , x\ = 1 U x 2 = . (63) 

The boundary of the phase space corresponding to the selection criterion 5 = 
7T — #ac is, in the limit m/ = 0, represented by 

x\-(l + x 2 )x l + — — X2 = {). (64) 

1 + cos 9 ac 

Therefore, if we require some acollinearity cut and M 2 (ff) > so, the limits of 
integration are: 

x 2 ml < x 2 < 1, x 2 < x x < 1, 

6 IFI interference contributions are generally small, as we know from pure O (a). However 
there is a recent analysis by S. Jadach and collaborators (private communication) which 
seems to indicate that exponentiation is quite important for the magnitude of interference 
contribution in cross-sections and asymmetries. We have been told that the difference between 
pure O (a) and exponentiation is quite often a factor 2. 



r,(i 



E+ = - (l-x 1 +x 2 ) v/i, E- = -xiy/a, E J = -(l-x 2 )Vs. (62) 



' ° <r rp <- ~acol 

S 



x 2 < £1 < x : 



*i" < *i < 1. 



(65) 



with 



acol 
•' 2 



1 - sin(6> aco i/2) 
1 + sin(0 aco i/2) : 
1 



l + x 2 ± 
2 



2 (l-2p)x 2 + l 



1 + cos 6> a 



1)1 



(66) 



The cut M 2 > 0.8 s deserves a comment. For a fully extrapolated setup 
we know from a complete calculation pj| that the correct scale in the coupling 
is s, i.e., a(s). Therefore the two codes adopt the following scales: loose M 2 
and s' cuts (up to 0.1s) with a(s), tight cuts with a(0). Now, however, we are 
using a tight M 2 -cut and there is an uncertainty, if one uses a(s) then partial 
higher corrections are seen, as expected. From Tab.(|35|) we see, however, that 
the net effect of IFI remains unchanged. The uncertainty associated with the 
choice of the scale, which becomes relevant only for high values of s'(M 2 )-cuts, 
is common to all final states, including leptons. 

Note that the T/Z agreement slightly deteriorates at such a tight cut (up to 
0.3 per- mill at resonance). This is not at all surprising: for such tight cuts some 
common exponentiation of ISR and FSR is mandatory. On the contrary, one 
should conclude that the T/Z agreement remains remarkable even for very tight 
M 2 cuts. 

In Tab.(p6[) we show the effect of changing the scale a{s) — > a(0) in the 
muonic cross-section and asymmetry with M 2 -cuts of 0.1 s, 0.5 s and 0.8 s. The 
variations in cross-sections range from rs 0.03 per-mill at M 2 = 0.1s up to 
« 4 per-mill at M 2 — 0.8 s. The effect of a scale change on the asymmetry is 
negligible. 

For the hadronic cross-section at M 2 = 0.8 the variation induced by a(s) — > 
a(0) is, at the peak, a 0.65 per-mill to be compared with 3.75 per-mill for muons; 
the sizeable difference is due to mass effects in QED FSR, m M <C m q . 
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Centre-of-mass energy in GeV 


M z -3 


M z - 1.8 


M z 


M z + 1.8 


M z +3 


No ISR/FSR interference, s' > 0.01 s 


o> [nb] 


0.22849 
0.22843 


0.47657 
0.47653 


1.48010 
1.47995 


0.69512 
0.69509 


0.40642 
0.40638 


+0.26 


+0.08 


+0.10 


+0.04 


+0.10 


CThad [nb] 


4.45012 
4.45146 


9.59910 
9.60165 


30.43639 
30.43824 


14.18269 
14.18391 


8.19892 
8.19923 


-0.30 


-0.27 


-0.06 


-0.09 


-0.04 


A» 

FB 


-0.28312 
-0.28330 


-0.16977 
-0.16985 


-0.00062 
-0.00066 


0.11186 
0.11182 


0.15466 
0.15461 


+0.18 


+0.08 


+0.04 


+0.04 


+0.05 


No ISR/FSR interference, s' > M 2 {^+^-) > 0.01 s 


o> [nb] 


0.22840 
0.22836 


0.47642 
0.47641 


1.47967 
1.47962 


0.69490 
0.69492 


0.40628 
0.40627 


+0.18 


+0.02 


+0.03 


-0.03 


+0.02 


Thad [nb] 


4.44990 
4.45129 


9.59865 
9.60132 


30.43501 
30.43725 


14.18203 
14.18342 


8.19853 
8.19894 


-0.31 


-0.29 


-0.07 


-0.10 


-0.05 


A£ 

FB 


-0.28321 
-0.28336 


-0.16981 
-0.16988 


-0.00062 
-0.00066 


0.11189 
0.11184 


0.15470 
0.15464 


+0.15 


+0.07 


+0.04 


+0.05 


+0.06 


ISR/FSR interference, s' > M 2 (fi+^-) > 0.01 s 


o> [nb] 


0.22783 
0.22779 


0.47566 
0.47566 


1.47970 
1.47967 


0.69563 
0.69570 


0.40684 
0.40686 


+0.18 


0.0 


+0.02 


-0.10 


-0.05 


CThad [nb] 


4.45083 
4.45221 


9.59990 
9.60254 


30.43495 
30.43717 


14.18082 
14.18215 


8.19761 
8.19797 


-0.31 


-0.27 


-0.07 


-0.09 


-0.04 


A» 

FB 


-0.28242 

-0.28287 


-0.16898 
-0.16934 


-0.00031 
-0.00031 


0.11177 
0.11193 


0.15467 
0.15473 


+0.45 


+0.36 


+0.00 


-0.16 


-0.06 



Table 32: Effects of ISR/FSR QED interference for cr M , cr had and A^ (in CA3- 
mode). First row T0PAZ0, second row ZFITTER third row relative (absolute) 
deviations in per-mill. The scale s is used for the running a. 
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Centre-of-mass energy 


in GeV 






M z -3 


M z - 1.8 


M z 


M z + 1.8 


M z +3 


u M [nb] 


T 


0.22311 


0.47122 


1.47473 


0.68963 


0.40076 






0.22153 


0.46914 


1.47464 


0.69138 


0.40204 






-7.13 


-4.43 


-0.06 


+2.53 


+3.18 


Z 


0.22315 


0.47127 


1.47466 


0.68967 


0.40079 






0.22157 


0.46922 


1.47463 


0.69148 


0.40211 






-7.13 


-4.37 


-0.02 


+2.62 


+3.28 


AL 


T 


-0.28347 


-0.16849 


0.00050 


0.11539 


0.16169 






-0.28074 


-0.16613 


0.00103 


0.11425 


0.16018 






+2.73 


+2.36 


+0.53 


-1.14 


-1.51 


Z 


-0.28344 


-0.16849 


0.00047 


0.11535 


0.16166 






-0.28084 


-0.16634 


0.00104 


0.11440 


0.16025 






+2.60 


+2.15 


+0.57 


-1.15 


-1.41 



Table 33: TOPAZO/ZFITTER comparison of IFI effect for M 2 > 0.8 s. First row 
is without IFI, second row is with IFI and third row is the net effect in per-mill. 
The scale is used for the running a. 





Centre-of-mass energy in GeV 




M z -3 


M z - 1.8 


M z 


M z + 1.8 


M z +3 


e>> [nb] 


T 


0.21932 

0.21775 

-7.21 


0.46287 
0.46081 

-4.47 


1.44794 

1.44785 

-0.06 


0.67725 
0.67896 

+2.52 


0.39366 

0.39492 

+3.19 


FB 


T 


-0.28450 
-0.28163 

+2.87 


-0.16914 

-0.16668 

+2.46 


0.00033 

0.00087 

+0.54 


0.11512 

0.11386 

-1.26 


0.16106 

0.15937 

-1.69 



Table 34: T0PAZ0 evaluation of IFI effect for M 2 (^+/i~) > 13.83997 GeV and 
#acoi < 10°. First row is without IFI, second row is with IFI and third row is 
the net effect in per-mill. The scale is used for the running a. 
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Chad [nb] 






Centre-of 


-mass energy in GeV 




M z -3 


M z - 1.8 


M z 


M z + 1.8 


M z +3 


T 


4.44990 


9.59865 


30.43501 


14.18203 


8.19853 


M 2 > 0.01s 


4.45083 


9.59990 


30.43495 


14.18082 


8.19761 




+0.21 


+0.13 


-0.002 


-0.09 


-0.11 


Z 


4.45129 


9.60132 


30.43725 


14.18342 


8.19894 


M 2 > 0.01s 


4.45221 


9.60254 


30.43717 


14.18215 


8.19797 




+0.21 


+0.13 


-0.003 


-0.09 


-0.12 


T 


4.36901 


9.45847 


30.05724 


13.97471 


8.05056 


M 2 > 0.8 s 


4.37172 


9.46205 


30.05739 


13.97167 


8.04831 


a(s) 


+0.62 


+0.38 


+0.005 


-0.22 


-0.28 


Z 


4.36943 


9.45848 


30.04915 


13.97099 


8.04804 


M 2 > 0.8 s 


4.37214 


9.46201 


30.04923 


13.96786 


8.04573 


a(s) 


+0.62 


+0.37 


+0.003 


-0.22 


-0.29 


T 


4.37186 


9.46463 


30.07686 


13.98432 


8.05646 


M 2 > 0.8 s 


4.37458 


9.46821 


30.07701 


13.98127 


8.05421 


a(0) 


+0.62 


+0.38 


+0.005 


-0.22 


-0.28 



Table 35: TOPAZO/ZFITTER comparison of IFI inclusion for hadrons with M 2 - 
cut. First entry is without IFI, second entry is with IFI, third entry is the IFI 
effect in per-mill. For the 0.8 cut there are two entries, corresponding to two 
different scales: a(s) and a(0). 
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Centre-of-mass energy in GeV 


M z -3 


M z - 1.8 


M z 


M z + 1.8 


M z +3 


o> [nb] 


T M 2 > 0.1 sa(s) 

T M 2 >0.1sa(0) 

Diff. 


0.22615 

0.22616 

0.04 


0.47369 

0.47370 

0.02 


1.47484 

1.47489 

0.03 


0.69179 

0.69182 

0.04 


0.40382 

0.40384 

0.05 


T M 2 > 0.5 sa(s) 

TM 2 > 0.5sa(0) 

Diff. 


0.22026 

0.22051 

1.13 


0.46327 

0.46380 

1.14 


1.44595 

1.44762 

1.08 


0.67712 

0.67792 

1.18 


0.39442 

0.39490 

1.22 


T M 2 > 0.8 sa(s) 

TM 2 > 0.8sa(0) 

Diff. 


0.20926 

0.21005 

3.76 


0.44215 

0.44381 

3.74 


1.38382 

1.38903 

3.75 


0.64513 

0.64769 

3.95 


0.37335 

0.37493 

4.21 


FB 


T M 2 > O.lsa(s) 

TM 2 > O.lsa(O) 

Diff. 


-0.28503 

-0.28505 

0.02 


-0.16986 

-0.16987 

0.01 


-0.00029 

-0.00029 




0.11224 

0.11225 

-0.01 


0.15557 

0.15558 

-0.01 


T M 2 > 0.5 sa(s) 

TM 2 > 0.5sa(0) 

Diff. 


-0.28603 

-0.28605 

0.02 


-0.16981 

-0.16982 

0.01 


0.00001 

0.00001 




0.11312 

0.11311 

0.01 


0.15727 

0.15726 

0.01 


T M 2 > 0.8sa(s) 

TM 2 >0.8sct(0) 

Diff. 


-0.28086 

-0.28091 

0.05 


-0.16613 

-0.16616 

0.03 


0.00130 

0.00128 

0.02 


0.11534 

0.11530 

0.04 


0.16236 

0.16225 

0.09 



Table 36: errand A^ B , including IFI, for various M 2 -cuts. First entry is with 
scale = s, second entry is with scale = 0, third entry is the relative (absolute) 
difference in per- mill. 



61 



7.2 Comparison with Realistic Kinematical Cuts 

Now we move to the classic setup, superimposing an acollinearity cut, a cut 
of 1 GcV on the energy of both outgoing fermions and a cut on the polar 
angular range of the outgoing fermion. The photon phase space is a little more 
complicated but not much. An acollinearity cut of 10° or 25° corresponds to 
M 2 = 0.84 s and M 2 = 0.64 s, and is thus regarded as a tight cut. Therefore, 
the scale for a is a(0). The comparison is shown in Tabs.(p7| Eof). 

There are several points to be discussed here. Although for the peak energy, 
we register generally a good agreement between TOPAZO 4 . 4 and ZFITTER 5 . 20, 
the situation for the off-peak points is signalling a disagreement. Here TOPAZO 
4 . 4 confirms the same size of IFI effect of previous situations, while ZFITTER 
5 . 20 sees half of the effect. This is a somewhat unique situation, the only 
one where TOPAZO 4.4 and ZFITTER 5.20 register a substantial disagreement. 
This disagreement should not be taken as a measure of the real theoretical 
uncertainty. 

Recent work of the Zeuthen group p0[ has contributed substantially in un- 
derstanding the origin of IFI-discrepancy, the IFI effect being so different in 
the two codes when realistic cuts are imposed. For the case of full angular 
acceptance, the work of |40) contains a full list of updated results. It recalcu- 
lates photonic corrections with acollinearity cuts, having in mind applications 
to ZFITTER. The conclusions of this work are that after inclusion of the new cal- 
culation the agreement with TOPAZO 4.4 is 0.1 per-mill (at the wings) or better 
(on resonance) for cross-sections. Work is in progress to update the inclusion of 
angular cuts. Thanks to the important contribution of the Zeuthen group, the 
IFI discrepancy with realistic cuts is on its way to be fully solved. 

As stated above, the introduction of angular cuts will represent the next step 
in Zeuthcn's program. For the moment, therefore, we have at our disposal only 
the comparison between TOPAZO 4.4 and ZFITTER 5.20. This section contains 
many IFI numbers based on these two programs, so that it should become clear 
under which conditions the discrepancies are large or small. Based on this the 
experiments will then have to derive an adequate solution for ROs with realistic 
cuts, concerning the magnitude of the effect and its theoretical uncertainty. 
Some possibilities are: 

1. Use current IFI as implemented in TOPAZO 4 . 4/ZFITTER 5 . 20 and assign 
overall uncertainty from the comparison, 

2. Use O (a) IFI of, e.g., K0RALZ to remove IFI effects from the ROs and fit 
the ROs with T0PAZ0/ZFITTER having switched off IFI. Of course, in this 
case the K0RALZ implementation of IFI must be evaluated to obtain the 
theory uncertainty. 

3. Wait for an updated ZFITTER code and repeat the T0PAZ0/ZFITTER com- 
parison to estimate the theoretical uncertainty. 
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7.3 Experimental Aspects of Initial-Final QED Interfer- 
ence 

As a last remark we observe that also the treatment of initial-final QED inter- 
ference effects by the experiments in arriving at quoted ROs has to be taken 
into account. For example, if the Monte Carlo generators used to correct for 
efficiency and to extrapolate for geometrical acceptance do not contain inter- 
ference (as is the case, for example, in the current versions of JETSET/PYTHIA, 
or KORALZ used in multi-photon mode), then the extrapolated and quoted ROs 
somehow miss interference. 

One main question is therefore: do QED interference effects enter only in 
the extrapolation step, e.g., cos(9) — ► 1? If so, then at least the fits are under 
control since the effect of missing IFI cancels out when TOPAZO/ZFITTER arc 
also run without interference. However, if interference effects already show up 
within the accepted region, then not only the quoted ROs are ill-defined but 
there are also problems in fitting them. 

It is of course a matter of the size of the effect. From our study we can 
draw some conclusion, despite the present disagreement in case realistic cuts 
are imposed. If we stay at the Z-peak where our predictions show agreement 
we observe a -0.07, -0.30, -0.68 (+0.01,-0.24, -0.65) pcr-mill effect in T0PAZ0 
(ZFITTER) on the cross-section for 6> acol < 10° and 8 acc = 0°,20° and 40°. The 
effect of interference grows when reducing the angular acceptance from 0° to 
40°. At the wings the effect is becoming smaller for reduced angular acceptance, 
but the rate of decreasing is not fast; for TOPAZO/ZFITTER the effect of IFI is still 
-3.67/ - 3.61 (1.18/1.22) per-mill at 6> acc = 40° and </s = M z - 3 (+3). Note 
that around 6> acc = 40° T0PAZ0 and ZFITTER start to agree well, so the previous 
statement is not affected by differences in the codes. For 6* aco i < 25° the effect 
at the peak goes from zero at 6> acc = 0° to —0.63 per-mill at 6* acc = 40°. At 
the wings the effect is —2 per-mill and +0.7 per-mill respectively for acc = 40° 
(with good agreement between the codes). For asymmetries, the interference 
effect decreases in magnitude with reduced angular acceptance at and below the 
peak energy. This is also the case for T0PAZ0 above the peak, while a different 
behaviour is observed for ZFITTER. 

Our conclusion is that the size of the IFI effect is at the level of few per-mill 
even within a fiducial volume cut. In order to minimise the effect we recommend 
to quote results, e.g., asymmetries, within a fiducial volume, at least at the wings 
of the Z resonance. 
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a M [nb] with aco i < 10° 






Centre-of- 


mass ener 


gy in GeV 




face 


M z -3 


M z - 1.8 


M z 


M z + 1.8 


M z +3 


T 0° 


0.21932 


0.46287 


1.44795 


0.67725 


0.39366 




0.21776 


0.46083 


1.44785 


0.67894 


0.39491 




-7.16 


-4.43 


-0.07 


+2.49 


+3.17 


z o° 


0.21928 


0.46285 


1.44781 


0.67722 


0.39361 




0.21852 


0.46186 


1.44782 


0.67814 


0.39429 




-3.48 


-2.14 


+0.01 


+1.36 


+1.72 


T 20° 


0.19990 


0.42207 


1.32066 


0.61759 


0.35886 




0.19873 


0.42049 


1.32027 


0.61873 


0.35973 




-5.89 


-3.76 


-0.30 


+1.84 


+2.42 


Z 20° 


0.19987 


0.42205 


1.32053 


0.61756 


0.35881 




0.19892 


0.42075 


1.32021 


0.61857 


0.35959 




-4.78 


-3.09 


-0.24 


+1.63 


+2.17 


T40° 


0.15034 


0.31762 


0.99428 


0.46479 


0.26989 




0.14979 


0.31680 


0.99360 


0.46514 


0.27021 




-3.67 


-2.59 


-0.68 


+0.75 


+1.18 


Z 40° 


0.15032 


0.31760 


0.99415 


0.46474 


0.26983 




0.14978 


0.31680 


0.99350 


0.46511 


0.27016 




-3.61 


-2.53 


-0.65 


+0.80 


+1.22 



Table 37: TOPAZO/ZFITTER comparison for the muonic cross-section of complete 
RO (CA3-mode) with the angular acceptance (# acc = 0, 20,40°) and acollinear- 
ity (^acoi < 10°) cuts. First row is without IFI, second row with IFI, third row 
is the relative (per-mill) effect of IFI. 
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a^ b with e acol < 10 


) 








Centre-of- 


mass energy in GeV 




"ace 


M z -3 


M z - 1.8 


M z 


M z + 1.8 


M z +3 


T 0° 


-0.28450 


-0.16914 


0.00033 


0.11512 


0.16107 




-0.28158 


-0.16665 


0.00088 


0.11385 


0.15936 




+2.92 


+2.49 


+0.55 


-1.27 


-1.71 


z o° 


-0.28453 


-0.16911 


0.00025 


0.11486 


0.16071 




-0.28282 


-0.16783 


0.00070 


0.11475 


0.16059 




+ 1.71 


+ 1.28 


+0.45 


-0.11 


-0.12 


T20° 


-0.27509 


-0.16352 


0.00042 


0.11171 


0.15645 




-0.27259 


-0.16146 


0.00084 


0.11064 


0.15499 




+2.50 


+2.06 


+0.42 


-1.07 


-1.46 


Z 20° 


-0.27506 


-0.16347 


0.00035 


0.11148 


0.15616 




-0.27408 


-0.16261 


0.00070 


0.11133 


0.15594 




+0.98 


+0.86 


+0.35 


-0.15 


-0.22 


T40° 


-0.24219 


-0.14396 


0.00054 


0.09906 


0.13903 




-0.24041 


-0.14262 


0.00078 


0.09837 


0.13809 




+ 1.78 


+ 1.34 


+0.24 


-0.69 


-0.94 


Z40° 


-0.24207 


-0.14386 


0.00050 


0.09893 


0.13891 




-0.24151 


-0.14343 


0.00069 


0.09890 


0.13888 




+0.56 


+0.43 


+0.19 


-0.03 


-0.03 



Table 38: T0PAZ0/ZFITTER comparison for the muonic forward-backward 
asymmetry of complete RO (CA3-mode) with the angular acceptance (0 acc = 
0,20,40°) and acollincarity (# aco i < 10°) cuts. First row is without IFI, second 
row with IFI, third row is the absolute effect of IFI (per- mill). 
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a M [nb] with # aco i < 25° 






Centre-of- 


mass energy in GeV 




"ace 


M z -3 


M z - 1.8 


M z 


M z + 1.8 


M z +3 


T 0° 


0.22333 


0.46971 


1.46611 


0.68690 


0.40034 




0.22233 


0.46838 


1.46611 


0.68812 


0.40127 




-4.50 


-2.84 


0.00 


+1.77 


+2.32 


z o° 


0.22328 


0.46968 


1.46598 


0.68688 


0.40031 




0.22281 


0.46905 


1.46603 


0.68754 


0.40081 




-2.11 


-1.34 


+0.03 


+0.96 


+1.25 


T 20° 


0.20359 


0.42835 


1.33731 


0.62648 


0.36507 




0.20286 


0.42733 


1.33700 


0.62725 


0.36569 




-3.60 


-2.39 


-0.23 


+1.23 


+1.70 


Z 20° 


0.20357 


0.42833 


1.33718 


0.62647 


0.36505 




0.20321 


0.42781 


1.33689 


0.62684 


0.36536 




-1.77 


-1.22 


-0.22 


+0.59 


+0.85 


T40° 


0.15320 


0.32245 


1.00698 


0.47167 


0.27479 




0.15286 


0.32190 


1.00635 


0.47185 


0.27500 




-2.22 


-1.71 


-0.63 


+0.38 


+0.76 


Z 40° 


0.15318 


0.32243 


1.00682 


0.47164 


0.27477 




0.15287 


0.32192 


1.00619 


0.47180 


0.27496 




-2.03 


-1.58 


-0.63 


+0.34 


+0.69 



Table 39: TOPAZO/ZFITTER comparison for the muonic cross-section of complete 
RO (CA3-mode) with the angular acceptance (# acc = 0, 20,40°) and acollinear- 
ity (^acoi < 25°) cuts. First row is without IFI, second row with IFI, third row 
is the relative (per-mill) effect of IFI. 
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A£ B with 9 acol < 25° 






Centre-of- 


mass energy in GeV 




Pace 


M z -3 


M z - 1.8 


M z 


M z + 1.8 


M z +3 


T 0° 


-0.28617 


-0.17037 


-0.00032 


0.11324 


0.15730 




-0.28501 


-0.16923 


0.00006 


0.11293 


0.15703 




+1.16 


+1.14 


+0.38 


-0.31 


-0.27 


z o° 


-0.28647 


-0.17049 


-0.00043 


0.11293 


0.15682 




-0.28555 


-0.16975 


-0.00005 


0.11307 


0.15701 




+0.92 


+0.74 


+0.48 


+0.14 


+0.19 


T20° 


-0.27695 


-0.16485 


-0.00026 


0.10974 


0.15250 




-0.27611 


-0.16404 


0.00001 


0.10955 


0.15236 




+0.84 


+0.81 


+0.27 


-0.19 


-0.14 


Z 20° 


-0.27722 


-0.16497 


-0.00037 


0.10944 


0.15204 




-0.27657 


-0.16447 


-0.00009 


0.10963 


0.15229 




+0.65 


+0.50 


+0.28 


+0.19 


+0.25 


T40° 


-0.24423 


-0.14536 


-0.00016 


0.09703 


0.13492 




-0.24381 


-0.14501 


-0.00004 


0.09702 


0.13500 




+0.42 


+0.35 


+0.12 


-0.01 


+0.08 


Z40° 


-0.24445 


-0.14545 


-0.00026 


0.09678 


0.13454 




-0.24444 


-0.14542 


-0.00011 


0.09700 


0.13483 




+0.01 


+0.03 


+0.15 


+0.22 


+0.29 



Table 40: TOPAZO/ZFITTER comparison for the muonic forward-backward 
asymmetry of complete RO (CA3-mode) with the angular acceptance (9 acc = 
0,20,40°) and acollincarity (# aco i < 25°) cuts. First row is without IFI, second 
row with IFI, third row is the absolute effect of IFI (per- mill). 
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8 Realistic Observables in the Model Indepen- 
dent Approach 

As we have shown in Eq.(|40]) the main emphasis in the Mi-approach is to or- 
ganise the calculation of ROs in terms of POs for MI fits. One has to show 
the internal consistency of the procedure. First, one has to show that for POs 
with values as calculated in the SM, the ROs are by construction identical to 
the full SM RO calculation. Next one has to make sure that the ROs computed 
by different codes also agree for MI calculations when the POs do not have SM 
values, but are varied over a range of PO values corresponding to at least the 
current experimental errors on POs for a single experiment. 

It is one of the main goals of our study to perform such comparison with 
TDPAZO and ZFITTER. Of course it will not be possible to present here the full 
comparison for arbitrary ROs but we will compare few significant examples, 
e.g., few ROs like cr M ,tr/j and A^ B , at fixed M z ,m t ,M H , a(M^),a s (M^), as a 
function of T z , cr®, Ri and A ' 1 . Note that we do not vary M z , both because its 
experimental error is so small and because it is also a SM parameter. We vary 
the remaining POs, one by one. Assuming the experimental errors to be twice 
the current LEP errors fujl, we obtain: 

T z = 2.4958 ±0.0048 GeV, 

a a h = 41.473 ± 0.116 nb, 

Ri = 20.748 ±0.052, 

A° F ' B = 0.01613 ±0.00192. (67) 

We vary the POs by twice their experimental errors. The corresponding devi- 
ations between the two codes will give an estimate of the technical precision 
of the procedure. For the SM initialisation, our preferred setup of SM input 
parameters is used. The effects of IFI and ISPP are not included. 

We show the MI comparisons in Figs. js|-|l5|. The solid curve gives devia- 
tions corresponding to the current experimental value of the relative PO, the 
dotted (dashed) curve corresponds to PO ± 2 x the experimental errors listed 
in Eq.([37]). For some of the figures the three curves become indistinguishable, 
typically the hadronic cross section has a constant deviation between the two 
Mi-calculations if we vary Ri or A® J . 

From Fig. || to Fig. [ll] we use an extrapolated setup with an s'-cut of 0.01 s 
and report the relative (in per- mill) T0PAZ0-ZFITTER deviations for 0^,07, and 
the absolute deviations for A^ B . In Figs. [l2j-[l5| we show similar comparisons for 
the muonic ROs with realistic cuts. We choose 20° < #_ < 160°, # a coii < 10°, 
and E t h (/x ) = 1 GeV. There is no appreciable difference in behaviour with 
respect to the s' case shown before. 
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Figure 8: Relative deviations between TDPAZO and ZFITTER for muonic and 
hadronic cross-section and absolute deviations for muonic forward-backward 
asymmetry in CA3 mode and for s' > 0.01 s. The solid curve gives deviations 
corresponding to the current experimental value of T z , the dotted (dashed) 
curve corresponds to PO ± 2 x the experimental error. 
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Figure 9: Relative deviations between TDPAZO and ZFITTER for muonic and 
hadronic cross-section and absolute deviations for muonic forward-backward 
asymmetry in CA3 mode and for s' > 0.01 s. The solid curve gives deviations 
corresponding to the current experimental value of <r[J ad , the dotted (dashed) 
curve corresponds to PO ± 2 x the experimental error. 
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Figure 10: Relative deviations between TDPAZO and ZFITTER for muonic and 
hadronic cross-section and absolute deviations for muonic forward-backward 
asymmetry in CA3 mode and for s' > 0.01 s. The solid curve gives deviations 
corresponding to the current experimental value of Ri , the dotted (dashed) curve 
corresponds to PO ± 2 x the experimental error. 
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Figure 11: Relative deviations between TDPAZO and ZFITTER for muonic and 
hadronic cross-section and absolute deviations for muonic forward-backward 
asymmetry in CA3 mode and for s' > 0.01 s. The solid curve gives deviations 
corresponding to the current experimental value of A®'^, the dotted (dashed) 
curve corresponds to PO ± 2 x the experimental error. 



72 



N 

I- 
< 

X 

CO 

O 



0.5 - 



- 



•0.5 - 



■1 



-i 1 r 



i — ■ — i — ' — i — ' — i — ' — i — ■ — r 



20° <9 < 160 c 



*aco. < 10 




A FB W 



_L 



J i I i L 



_L 



J i L 



87 88 89 90 91 92 93 94 95 



0.5 



CD 

N 

I- 



- 



•0.5 



■1 



T z variation 




plus 

minus 

central 



a(/i.) 



_L 



_L 



_L 



_L 



_L 



_L 



87 88 89 90 91 92 93 94 95 

E rm [GeV] 

cm 

Figure 12: Relative deviations between TOPAZO and ZFITTER for muonic cross- 
section and absolute deviations for muonic forward-backward asymmetry in CA3 
mode and for realistic cuts. The solid curve gives deviations corresponding to 
the current experimental value of T z , the dotted (dashed) curve corresponds to 
PO ± 2 x the experimental error. 
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Figure 13: Relative deviations between TOPAZO and ZFITTER for muonic cross- 
section and absolute deviations for muonic forward-backward asymmetry in CA3 
mode and for realistic cuts. The solid curve gives deviations corresponding to 
the current experimental value of c[j ad , the dotted (dashed) curve corresponds 
to PO ± 2 x the experimental error. 
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Figure 14: Relative deviations between TOPAZO and ZFITTER for muonic cross- 
section and absolute deviations for muonic forward-backward asymmetry in CA3 
mode and for realistic cuts. The solid curve gives deviations corresponding to 
the current experimental value of Ri , the dotted (dashed) curve corresponds to 
PO ± 2 x the experimental error. 
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Figure 15: Relative deviations between TOPAZO and ZFITTER for muonic cross- 
section and absolute deviations for muonic forward-backward asymmetry in CA3 
mode and for realistic cuts. The solid curve gives deviations corresponding to 
the current experimental value of A®^ , the dotted (dashed) curve corresponds 
to PO ± 2 x the experimental error. 
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Most of the plots in Figs. p|~[l5| show a crossing of the three curves (central, 
plus and minus) at s = M 2 z . This means that at the Z-pole the MI imple- 
mentations of the two code are fully equivalent while the MI treatment of the 
off-resonance terms is somewhat different. This fact is largely expected since 
Mi-implementations or the MI-SM splitting is far from unique. However, the 
rather satisfactory level of agreement is telling us that the associated theoretical 
uncertainty in extracting POs from ROs is not substantially different or badly 
deteriorated with respect to the one that we have shown in SM comparisons. 

Another typical effect is that TOPAZO ROs tend to be lower than ZFITTER ROs 
at PO — 2 x experimental error and on the low-energy side of the resonance. 
With increasing energy the curves tend to cross at the Z-peak and to reverse 
their sign on the high energy side. On the contrary, for PO + 2 x experimental 
error, TOPAZO is higher below the resonance and lower above it. 

In Fig. [l(| we compare TOPAZO and ZFITTER predictions by computing a^, ah 
and A^ B with and s'-cut of 0.01 s. The solid line gives deviations in the SM 
predictions initialised with our preferred setup. The dotted line gives deviations 
from the two MI predictions in case the ROs are computed in term of POs 
evaluated at their (code dependent) SM values. From comparing the two sets of 
curves we see that no serious degradation arises in the transition RO(SM,T,Z) 
— > RO(POgM, T, Z). In particular we continue to register a very good agreement 
around the peak. 

From these figures it emerges that, compared to the SM comparison, the 
agreement is still reasonable. Much more interesting than the deviations shown 
here would be to compare differences between MI fits performed with TOPAZO 
and ZFITTER using the same RO data set, clearly a task to be done by the 
experimental collaborations. The differences shown in the figures will cause 
differences in Mi-fit results between TOPAZO and ZFITTER for the same input set 
of ROs. However, the effects have to be seen compared to the fit (experimental) 
errors. Approximately, the differences in ROs seen at the pole centre-of-mass 
energy will be those observed in the corresponding POs. 
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Figure 16: Relative deviations between TDPAZO and ZFITTER for muonic and 
hadronic cross-section and absolute deviations for muonic forward-backward 
asymmetry in CA3 mode and for s' > 0.01 s. The solid curve gives deviations 
corresponding to SM predictions, the dotted curve corresponds to a MI predic- 
tion where PO arc initialised to their SM values. 
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9 Theoretical Uncertainties for Realistic Observ- 
ables 

When varying options in the calculation of ROs based on POs or SM parame- 
ters, the theoretical uncertainties on ROs are obtained. This is discussed in the 
following. However, when the experimental collaborations analyse their mea- 
surements of ROs, this procedure is inverted by fitting POs or SM parameters to 
the measured ROs: Varying options, the changes in fitted parameters reflect the 
theoretical uncertainties associated with the calculation of ROs. This way the 
theoretical uncertainties on ROs are propagated back to the fitted parameters. 
Care must be taken not to double-count theoretical uncertainties in the SM cal- 
culation of POs which propagate to ROs, when SM parameters are determined 
in a fit to POs themselves determined in a MI fit to ROs. 

9.1 Uncertainties in Standard Model Calculations 

We now propagate the various electroweak options in SM calculations fr om 
PO = PO(SM) to RO = RO(SM). These options are described in Section U 
and are those used in describing the theoretical uncertainties at the level of 
PO. In addition there are other uncertainties, which we have already discussed, 
like those associated to different treatments of initial-state QED radiation. The 
corresponding flag values in TDPAZO/ZFITTER are discussed in the following. 

In TOPAZO the flag 0HC='Y' selects next-to-leading and higher-order hard- 
photon contributions. Therefore this flag should not be varied. OnceOHC='Y' is 
initialised TOPAZO 's flag ORAD will select the type of next-to-leading and higher- 
order hard-photon contributions to be included, (A,D,E,F,Y). F J3J] is the 
recommended choice; Y ]3q | has been implemented after version 4.3 and uses 
the order a 3 YFS radiator. The other choices, A,D,E, remain for compatibility 
tests with previous versions and are not to be included in the estimate of the 
theoretical uncertainty. 0FS= ' D ' , ' Z ' selects the treatment of higher-order final- 
state QED corrections, see below. For the remaining TOPAZO's flags we observe 
the following. With OWBOX weak boxes are included ' Y ' or not ' N ' . The correct 
choice is always 'Y'. At the Z-resonance, but only there, weak boxes can be 
neglected having a relative contribution < 10 -4 . To give an example, for the 
muonic cross-section the effect of weak boxes is —0.01 per-mill at s/s = M z 
and +0.17 (-0.15) per-mill at the left (right) wing. Flags OAAS, DWEAK must 
be kept equal to ' N ' , ' R ' , respectively. 

In ZFITTER, the flag F0T2 controls which radiator is used for inital-state QED 
radiation, F0T2=-1 ,0,1,2,3,4,5. Values F0T2=3 , 5 switch between the additive 
(default) and the factorized order a 3 radiator, respectively. The remaining F0T2 
flag values are kept for compatibility tests with previous versions and are not 
to be included in the estimate of the theoretical uncertainty. Weak boxes are 
switched off/on with the flag B0XD=0,1, respectively, and must always be on. 
The flag CDNV controls the convolution of electroweak couplings and is discussed 
in Section g. This flag must be kept fixed at a value larger than zero. 
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For the muonic and hadronic cross-sections and for the muonic forward- 
backward asymmetry we find the results reported in Tab.(ETl). These results 
are obtained by running TOPAZO under different options. As for POs, we show 
the central value for ROs evaluated at the preferred setup; the minus error 
for RO central - min opt RO and the plus error for max opt RO - RO ce „ trai . For 
cross-section the errors are reported in pb and, when available, we also show 
the absolute differences T-Z. From Tab. (Ell) we observe that the difference be- 
tween TOPAZO and ZFITTER RO is fully compatible with the TOPAZO-estimated 
theoretical uncertainty. 

Finally we illustrate the effect of different treatments of final-state QED radi- 
ation in the presence of severe kinematical cuts. A possible source of theoretical 
uncertainty can be introduced when cuts are present, due to a different treat- 
ment of higher-order final-state QED effects: it can lead to differences which in 
general depend on the experimental cuts required and that may grow for par- 
ticularly severe cuts, ft was already shown jjj that two possible prescriptions, 

• completely factorized final-state QED correction versus 

• factorized leading-terms and non-leading contributions summed up, 

can lead to substantial differences. In TOPAZO the flag OFS selects the treat- 
ment of higher-order final-state QED corrections; (D) or (Z), respectively. In 
ZFITTER, the first option, completely factorized final-state QED correction, is 
implemented. TOPAZO predicts for er^ and 20° < (9_ < 120°, E th > 15GeV, 
#acoi < 10° an uncertainty of 1.2 per-mill at the five energy points. On the 
contrary for a loose cut of Eth > 1 GeV we obtain a reduction to a mere 0.1 
per-mill. 

9.2 Uncertainties in Model Independent Calculations 

Having discussed the theoretical uncertainties associated with the SM calcu- 
lations of ROs, we also need to address the question of uncertainties for the 
Mi-calculations, RO = RO(PO). The result is that one has all errors already 
quoted for RO(SM), Tab.(fll|), those derived from varying PO away from 
their SM values and obtained by comparing TOPAZO with ZFITTER, Figs. |)-^5[ 

Note that in the SM there are relations like POi = /(P0 2 ) = F(M Z , m t , . . .)• 
With PO = PO(SM) one has to make sure that, no matter how the Mi-structure 
is built, MI(T) must be equal to MI(Z). But when we break the SM-relations 
then several possibilities arise: one may write everywhere POi or /(PO2) and 
the difference is a measure of the associated uncertainty. 

Alternatively one can estimate this uncertainties internally to each code by 
running in MI mode with all the relevant electroweak and QED flag variations 
as discussed before. The structure of the MI calculations is such that, given 
the decomposition PO SM remnant, electroweak flag changing will effect the 
latter (the SM complement) and leave unaltered the PO (MI) component. To 
give an example, the whole Z — 7 interference in the quark sector is taken from 
the SM remnant and, therefore, influenced by flag setting. 
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0.04 pb 




^ 


M z -3 


0.14659 nb 


0.02 pb 


< 0.01 pb 




M z - 1.8 


0.27874 nb 


0.03 pb 


< 0.01 pb 




M z 


0.74051 nb 


0.04 pb 


0.04 pb 




M z + 1.8 


0.30868 nb 


0.02 pb 


< 0.01 pb 




M z +3 


0.17178 nb 


0.02 pb 


< 0.01 pb 




Chad 


M z -3 


4.45012 nb 


0.99 pb 


1.40 pb 


-1.29 pb 


M z - 1.8 


9.59909 nb 


1.81 pb 


3.41 pb 


-2.49 pb 


M z 


30.43639 nb 


1.85 pb 


14.27 pb 


-11.83 pb 


M z + 1.8 


14.18269 nb 


2.14 pb 


6.01 pb 


-1.27 pb 


M z +3 


8.19892 nb 


1.46 pb 


3.38 pb 


-0.36 pb 



Tabic 41: Theoretical uncertainties for G a . A^ 1 



and for cr had from T0PAZ0. 
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10 Production of Secondary Pairs 

Finally we come to the inclusion of pair- production in the calculation of realistic 
observables. Radiative photons from the initial- or final-state fermions may 
convert, leading to additional (soft) // pairs besides the primary pair. This 
leads to the problem of the signal definition, i.e., what is considered as (radiative 
correction to) fermion-pair production, and what is considered as genuine four- 
fermion production. 

Note that most Monte Carlo event generators used for fermion pair pro- 
duction do not include the radiative production of secondary pairs. In case of 
final-state pair production, visible in the detector, this may bias efficiency and 
acceptance calculations, in particular if the primary pair is a lepton pair. 

10.1 Initial-State Pair Production 

A fermionic pair of four-momentum q 2 radiated from the e + or e~ line gives 
a correction which is computed in [|4lfl . Also for this term there are different 
treatments, i.e., we can exponentiate the pair-production according to the YFS 
formalism |4^] or the same pairs can be included at O (a 2 ) . The physical uncer- 
tainty on the pair correction is given by that on the contribution of light-quark 
pairs, and is estimated to be 1.8 ■ 10~ 4 for cross-sections J42J. 

Both T0PAZ0 and ZFITTER adopt a hybrid solution which gives a remarkable 
agreement around the peak with the results of []42|. ZFITTER includes the radi- 
ation of e, /i, r-pairs and hadronic pairs. TDPAZO does not include radiation of 
r-pairs, their contribution is below the accuracy requirement of O (lO -4 ). For 
technical details we refer to since no further upgrading has been performed 
in the area of pair production. 

A cut was selected so that z m [ n s = s' = M 2 (ff) > 0.25 s, where // denotes 
the primary pair, i.e., z m i n is the minimum fraction of squared invariant mass 
of the final state (primary pair) after radiation of the additional initial-state 
pair (secondary pair). Since we neglect terms coming from ISPP (g> FSR there 
will be no difference between s'-cuts and M 2 -cuts as far as ISPP corrections are 
concerned. The soft-hard separator A has been fixed in the region where we see 
a plateau of stability. 

The effect of including leptonic as well as hadronic pairs radiated from the 
initial state is summarised in Tab.([42"|). From Tab. (|42]) we observe that the in- 
clusion of initial-state pair-production lowers the cross-sections up to an energy 
of M z + 3 where the effect becomes positive. At the peak pair production mod- 
ifies a^ by -2.55 (-2.53) per-mill according to T0PAZ0 (ZFITTER) and cr ha d by 
-2.55 per-mill for both codes. 

10.2 Final-State Pair Production 

The current versions of T0PAZ0 (4.4) and ZFITTER (5.20) do not include effects 
of final-state pair production. However, virtual and real corrections due to 
final-state pair production cancel to a large extent and the remaining effect 



N2 







Centre-of- 


mass energy in GeV 




M z - 3 


M z - 1.8 


M z 


M z + 1.8 


M z +3 


CT/i 


T 


0.22849 


0.47657 


1.48010 


0.69512 


0.40642 




0.22796 


0.47534 


1.47633 


0.69480 


0.40713 




-2.32 


-2.58 


-2.55 


-0.46 


1.75 


Z 


0.22843 


0.47653 


1.47995 


0.69509 


0.40638 




0.22790 


0.47532 


1.47621 


0.69478 


0.40708 




-2.33 


-2.55 


-2.53 


-0.45 


1.73 


FB 


T 


-0.28312 


-0.16977 


-0.00062 


0.11186 


0.15466 




-0.28377 


-0.17020 


-0.00062 


0.11192 


0.15439 




-0.65 


-0.43 


0.00 


+0.06 


-0.27 


Z 


-0.28330 


-0.16985 


-0.00066 


0.11182 


0.15461 




-0.28395 


-0.17028 


-0.00066 


0.11187 


0.15434 




-0.65 


-0.43 


0.00 


+0.05 


-0.27 


Chad 


T 


4.45012 


9.59910 


30.43639 


14.18269 


8.19892 




4.43937 


9.57410 


30.35866 


14.17587 


8.21336 




-2.42 


-2.60 


-2.55 


-0.48 


1.76 


Z 


4.45146 


9.60165 


30.43824 


14.18391 


8.19923 




4.44070 


9.57659 


30.36069 


14.17709 


8.21345 




-2.42 


-2.62 


-2.55 


-0.48 


1.73 



Table 42: The effect of including initial state pair production in T0PAZ0 and 
ZFITTERin CA3-mode, s' > 0.01 s. First (fourth) entry is without, second (fifth) 
entry is with pair production. Third (sixth) entry is the net effect in per-mill 
of the inclusion. 



is mostly absorbed in the running electromagnetic coupling a(s) entering the 
correction factor for FSR p3| . Therefore, the experimental event selections 
should not discriminate against additional (soft) pairs. Otherwise, a correction 
has to be applied before T0PAZ0/ZFITTER calculations can be compared with 
the measurements. 



n:-5 



11 Conclusions 

In Tab. (Ell) we show an estimate of the theoretical error for realistic observables 
as computed internally by TDPAZO. Another piece of information is given by the 
differences TOPAZO - ZFITTER among the (theoretical) central values for each 
quantity: these differences are basically (even though not totally) a measure of 
the effect induced by a variation in the renormalization scheme. 

Consider, however, the complete hadronic cross-section (in CA3 s'-mode): 
the differences (in per-mill) are -0.30, -0.27, -0.06, -0.09, -0.04 for the five 
centre-of-mass energies. There are two different origins for them, differences 
already present in the de-convoluted cross-sections (SD-mode) and differences 
due to convolution with initial-state QED radiation. For the former we find (in 
per-mill) 

-0.38, -0.30, -0.07, -0.15, -0.14. (68) 

The effect of convolution is found to be (in per-mill) 



0.06, +0.02, +0.01, +0.07, +0.14. 



(69) 



Especially on the high-energy side of the resonance we observe a partial com- 
pensation, leading to a small overall uncertainty. A more conservative attitude 
consists in adopting some rough approximation thus defining 



and adding the errors in quadrature: 



SD 



(1 + S dcc ) 



(70) 



Act 

dcc\2 



Act 



SD 



SD,T 



,SD,Z 



"Act" 
ct sd_ 


2 

+ 


FctAct sd i 

1(ctSD) 2 J 



AS 



In this way we end up with an estimate of the theoretical error of 

2.73 pb, 4.76 pb, 3.51 pb, 3.23 pb, 1.71 pb, or 
0.061%, 0.050%, 0.012%, 0.023%, 0.021%, 



(71) 



(72) 



for the complete hadronic cross-section. Although quite conservative, we con- 
sider the above as a safe estimate of the theoretical error. A comparison with 
the results of Tab. (El]) shows a substantial agreement with the estimate made 
internally by TOPAZO. For the complete muonic cross-section we find 

0.07 pb, 0.06 pb, 0.19 pb, 0.11 pb, 0.08 pb or 

0.030%, 0.014%, 0.013%, 0.016%, 0.021%. (73) 

In both cases, the uncertainty arising due to the uncertainty on the ISPP con- 
tribution has to be added. 

For the muonic forward-backward asymmetry we find results which are prac- 
tically indistinguishable from those already reported in Tab. (Ell). 
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